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ABSTRACT 


The cO., laser induced breakdown of Ar under the influence of an 
arc discharge is investigated. The laser beam is focused inside and 
outside the discharge. Inside the discharge the changing electron 
density and the reduced neutral density modify the breakdown thresholds. 
Outside the discharge, the effects of uv preionization and shock wave 
expansion on the breakdown thresholds are investigated. 


The pulsed CO, laser system, the electric discharge, and the break- 


2 


down chamber are described. Langmuir probes, interferometric and CO, 
absorption measurements are performed on the discharge plasma which 
yield information on the electron and neutral atom concentrations. 
Rotating Langmuir probes are used to measure uv photoion densities out- 
side the discharge. Shock wave strength measurements using piezo- 
electric probes are also reported. 

A cascade model of the breakdown process is developed. Numerical 
calculations of the plasma build up in the laser focus are reported. 

The experimental thresholds are compared with the results of the 
calculations. Good agreement is obtained. This leads to the possibi- 
lity of using the breakdown phenomenon as a diagnostic for electron 


and neutral atom components of a plasma. The feasibility of this idea 


is discussed. 
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CHAPTER I 


INTRODUCTION 


One of the most interesting aspects of the interaction of intense 
electromagnetic fields with matter is the breakdown of gases under the 
influence of a strong laser light. Gases which are normally transparent 
to light at low laser intensity become strongly ionized at some value of 
the intensity called the breakdown threshold ,resulting in strong absorpt- 
ion of the incident laser light and breakdown of the gas. The plasma 
produced expands rapidly and emits light and, hence, it is called a "spark". 
A brief review of the literature on laser induced gas breakdown with 


emphasis on CO, laser breakdown can be found in Chapter II. 


2 
Laser induced breakdown has attracted considerable attention in 
recent years. Numerous investigations have been performed both experiment- 

ally and theoretically which have provided much information on the nature 
of the phenomenon. The dependence of the breakdown threshold on gas 
pressure, focal length, ionization potential, pulse length and other para- 
meters have already been investigated. In the case of moderate to high 
pressures,the breakdown thresholds have been explained by the cascade 
theory of laser induced gas breakdown which says that the electon-ion 
density produced by a focused laser beam increases exponentially through 

an avalanche process, i.e., each electron present in the focal volume gains 
energy from the laser field via inverse bremsstrahlung and creates another 
electron through an ionizing collision with a gas atom or molecule. The 


two electrons then repeat the above process until a sufficiently high 


electron density is reached at the focus. 
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The cascade theory assumes the existence of initial or priming 
electrons at the focal volume to start the avalanche at the Start of the 
laser pulse. The more electrons there are at the focus, the fewer 
generations of electrons needed for breakdown to occur and hence, the smaller 
the required laser power. Breakdown measurements with deliberate pre- 
ionization of the gas have been performed and indeed, a lowering of the 
threshold has been observed. 

The primary motivation for doing this project is to determine the 
dependence of the breakdown threshold on the initial preionization and 
attempt to explain the results using cascade theory. Although cascade 
theory has adequately explained breakdown threshold measurements at high 
pressures, all experiments so far have only been "global" verifications 
of the theory and there are still some who doubt its validity. A success- 
ful cascade theory prediction of the dependence of the threshold on the 
initial electron ieustey would be strong evidence in its favor. To this 
end, a co, laser is focused in Ar which has been preionized by an electric 
discharge. The laser is focused not only in the region filled by the discharge 
but also in the uv preionized volume away from it. This provides a broad 
range of initial conditions needed in the study. The effect of temperature 
on the breakdown threshold can also be studied in this system. 

The dependence of the threshold on gas pressure or equivalently, the 
neutral density, is well known. Hence, if the electron density of a 
partially ionized gas is known, it may be possible to use the breakdown 
phenomenon tomeasure neutral gas density. Gas breakdown as a diagnostic 
tool would be very useful because it could provide a very high degree of 


space and time resolution, these being determined by the focal spot size and 
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the laser pulse length, respectively. At the same time, it would be a 
very convenient diagnostic tool because we would not need very sophisticated 
detection systems. This work will investigate the possibilities of gas 
breakdown being employed as a diagnostic tool. 

The exact eonditions in the focus at breakdown are not known. 
There can be a multitude of possible definitions that can and have been 
used in various experiments. These include the onset of light emission or 
laser absorption, visual observation of the spark and so on. Theoretically, 
threshold can be defined as the laser power which will give a certain 
value of the electron-ion density at the focus. This problem is discussed 
in greater detail in Chapter III. 

A description of the various experimental procedures and setups used 
are presented in Chapter IV. The various diagnostics performed on the uv 
photoplasma and the electric discharge are presented in Chapters V and VI. 
Chapter V also includes a discussion of probe theory which attempts to extend 
existing probe theories for high pressure flowing plasmas reported in the 
literature. In Chapter VII, a cascade model of gas breakdown is constructed 
and numerical calculations are performed which adequately explain the 
experimental results. The measured breakdown thresholds are reported in 
Chapter VIII. A comparison with cascade theory predictions are made. The 
possibility of using the breakdown phenomenon as a diagnostic for electron 
and neutral densities is explored. Finally, Chapter IX contains a discussion 


and summary of the important results of this investigation. 
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CHAPTER II 


REVIEW OF PREVIOUS LASER INDUCED GAS BREAKDOWN WORK 


2.1 Outline 

A brief review of previous work on laser induced gas breakdown 
is presented. Multiphoton and cascade ionization as possible electron 
production mechanisms occurring in the laser focus are discussed. The 
dependence of the breakdown threshold on pressure, focal length, pulse 
length, ionization potential, frequency, and other parameters as 
predicted by theory and as measured by various investigators are summarized. 
Cascade theory is seen to apply to gas breakdown at high pressures 


when using long laser pulses. 


2.2 Stages of Gas Breakdown 


The first observation of laser induced gas breakdown was reported by 
P.D. Maker, R.W. Terhune, and C.M. Savage in 19637. This was subsequently 
followed by vigorous research by many investigators. Reviews of laser 
induced breakdown work up to 1971 are available in the literature’ 

The breakdown phenomenon is usually divided into four stages as 
follows: 

1. preionization or initiation 

2. cascade stage 

3. plasma development and expansion 


4. extinction 


In this work, only the cascade stage will be discussed in any detail. The 
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other stages are thoroughly discussed in the review papers quoted above. 

A good discussion of the initial stages of breakdown is given by Papouler 
A striking feature of laser induced gas breakdown is the sudden 

change in the opacity of the gas to laser light. This is due to the 

sudden increase in the plasma density at the focal volume in the presence 

of the intense optical field. Figure 2.1 shows typical laser pulse 

shapes transmitted with and without breakdown. Two processes that can 

explain the ionization by high intensity light are multiphoton ionization 


and cascade ionization. 


Je Multiphoton Ionization 


In multiphoton ionization, several photons are absorbed simultaneously 
by a bound electron resulting in its detachment from the parent atom. 
This is the ionization mechanism believed to operate during the preioniza- 
tion stage of gas breakdown. At this time gas atoms or impurities with 
low ionization potentials are directly ionized by the laser electric field. 
The theory of multiphoton ionization has been developed into an 
elegant mathematical theory by numerous euthore! ius All multiphoton 
theories basically lead to a common prediction. The number of ions N, 
produced by light of intensity I is 


N 


ro) 
N, a No (Ihv) 


where N. is the neutral atom density and Ny is the number of photons whose 
a 
total energy is equal to or greater than the ionization potential of the 


gas. Experimentally, multiphoton ionization can be observed by going to 
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very low pressures so that the electron mean free path is greater than 
the focal spot diameter or by using ultrashort laser pulses such that the 
pulse length is shorter than the characteristic cascade time. In all 
cases, theoretical predictions are orders of magnitude higher than 
experimental results. This may be due to smearing of levels due to the 


very high local fields resulting in quasi-resonant transitions. 


2.4 Cascade Ionization 

In the case of high pressures or long laser pulses cascade ionization 
is believed to occur in the focus. The first electrons produced either 
by the multiphoton process or by deliberate preionization of the gas 
(preionization stage) can absorb energy from the laser field during 
collisions with neutral atoms via inverse bremsstrahlung. The electrons 
that acquire kinetic energies greater than the ionization potential of 
the gas can ionize the gas during collisions with gas atoms. The electron 


density increases exponentially: 


where y = Qn2 is the cascade development constant and t is the time 


al|e 


it takes to double the electron concentration. A more detailed discussion 
of cascade theory can be found in Chapter VII. References 12-15 also 
give good discussions of cascade theory and its basic results. For now 


we summarize its important features. 


i The rate the electrons gain energy from the field is the 


same as that in microwave theory, namely 
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dt m 2 Z 
+ 
sete 
Here Vofe = LPs is the collision frequency, n. is the 


neutral atom density, v is the electron velocity, oe is 

the collision cross-section for momentum transfer, e and m 

are the charge and mass of the electron and E and w are 

the intensity and frequency of the laser electric field. 

The significant loss mechanisms during the electron build-up 
are electron diffusion out of the focal volume, elastic and 
inelastic collisions, and at high pressures, electronic re- 
combination. Excitation may sometimes be a gain process 
depending on the closeness of the excited energy level to 

the ionization continuum, the laser intensity and the photon 
energy or when the excited level is a quasi-resonant level. 
Diffusion loss dominates at very low pressures and is manifested 
as a characteristic ae dependence of the threshold on 
pressure. At moderate pressures and long pulses, the threshold 
is determined by the laser pulse length and the initial elect- 
ron density and is seen as a pee dependence of the threshold 
on gas pressure. The elastic loss term is independent of 
pressure and becomes important only at very high pressures. 
Excitation loss is more important at lower laser fluxes and 
since the threshold decreases with increasing pressure, its 
effect is relatively greater at high pressure. Excitation thus 


weakens the pressure dependence of the threshold especially 
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in the case of gases with many excited levels such as 


molecular pases 13,18-19 


Numerical calculations taking 
into account these gain and loss mechanisms have been 
performed and have successfully accounted for measured break- 


down thresholds. 


2.5 Experimental Results 
Dependence of the breakdown threshold on gas pressure 


The pressure dependence of the threshold has been determined using 


ruby, neodymium and CO, decere uno Attention has been focused on noble 


gases, nitrogen, and CO For example, Meyerand and heen. studied 


2s 
breakdown of He and Ar. In Figure 2-2 the results of Hill, James and 


Peredena: using a CO, laser are reproduced showing the typical pressure 


2 
dependence of the threshold. The slopes in the decreasing portion of the 
curves are between -l and -2 as predicted by cascade theory. The 
existence of the threshold minimum when w = Vofe is also predicted by 
cascade theory. A high magnetic field will tend to reduce the diffusion 
loss out of the focal volume. Hence it is expected that there will be a 


decrease in the breakdown threshold upon application of an external magnetic 


LS 
field. This was observed by Cohn, et.al... 


Dependence of the threshold on the ionization potential of the gas 


Experimental results Bhowaa lee that the threshold generally increases 


with increasing ionization potential. This is predicted by both cascade 
and multiphoton theories. The only exception is Ne which has a higher 


breakdown threshold than He. 
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Dependence of the breakdown threshold on the radiation frequency 

The frequency dependence of the breakdown threshold has been studied 
by Buscher, et. Tikes Akhmanov, et. ine righ oe and Morgan, et. aiibee 
Buscher, et. al. and Akhmanov, et. al. used 1.06 um radiation from a 
Nd glass laser and its harmonics. Morgan, et. al. used ruby, Nd, and 
a tunable dye laser. It is generally observed that the threshold increases 
with frequency as ie in agreement with theory. However, at still higher 
frequencies, the threshold starts to decrease. This might be due to 


enhancement of the photoionization of the excited atoms and an increase in 


the probability of producing the priming electrons from impurities. 


Dependence of the threshold on the mode pattern of the laser output 


A multimode laser when focused will produce a very inhomogeneous 
intensity distribution. Some areas in the focal spot will have intensities 
up to 2 orders of magnitude higher than the mean intensity. For a Gaussian 
mode pattern, on the other hand, the maximum value of the intensity at 
the center of the spot is only 1.4 times the average value. Smith and 
poninson-© and Alcock, et. Ape have found no significant difference in 
the breakdown threshold when one uses a single mode or multimode laser. 
This suggests that the diffusion processes are effective in averaging out 


the effect of the inhomogeneities and breakdown occurs under the effect of 


the average field. 
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Dependence of the breakdown threshold on the focal volume 


The focal length dependence of the threshold has been studied and 
reported in references 28-30. The threshold is a function of the focal 
length since the focal spot diameter and therefore the diffusion length 
is proportional to it. However, even at high pressures, the threshold 
has been observed to depend on the focal length. Experimental results 
also tend not to agree with theoretically predicted thresholds. It is 
believed that this anomalous focal spot size dependence is due to the 
presence of dust or other impurities which are easily ionized. The 
larger the focal spot, the greater the probability of finding an impurity 


particle in it and hence, the lower the breakdown threshold. 


Dependence of the breakdown threshold on the laser pulse length 

The dependence of the threshold on the pulse length is determined 
by the dominant loss mechanism. For short laser pulses, diffusion loss 
can be neglected and the threshold varies inversely as the pulse length. 
For very long laser pulses, the: threshold is practically independent of 
the pulse length, depending mainly on the peak intensity. The pulse length 
dependence of the breakdown threshold has been experimentally determined 
29,31, 32 


by several investigators for long and short pulses Measurements 


by Hacker, Cohn and tare point to the conclusion that cascade theory can 


describe the breakdown process. 


Dependence of the breakdown threshold on the initial electron density 


Little work has been done on this aspect of laser induced gas break- 


down to the present. Breakdown thresholds in uv preionized Ar were 
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performed but the thresholds were not related to the actual initial 


5... 39=34 
electron concentration - Others have focused the laser right into a 


355,30 


discharge In this case the initial electron concentration can be 


estimated from the current and voltage across the discharge if the electron 
temperature is assumed to be equal to the gas temperature which was taken 
to be 300 Ree In all measurements with preionization, a depression of 

the breakdown threshold was observed with increasing electron preionization. 
This was attributed to a transition from free to ambipolar diffusion loss 
during the cascade. Also, the higher the initial electron density, the 
lower number of electron generations required before breakdown as would 

be expected for a cascade process. A closer investigation of this aspect 
of laser induced breakdown and in particular breakdown of uv preionized 


Ar is the subject of the present work. 


2.6 Expansion and Extinction Stages 


After the occurrence of breakdown, the plasma becomes highly absorbing 
and is rapidly heated up by the incident laser pulses. Hydrodynamic 
expansion and shock expansion now take place. The spark tends to move toward 
the direction of the laser aside from the usual radial expansion. The 
expansion towards the laser is similar to a detonation wave. At the end of 
the pulse, the plasma cools down and decays mainly by recombination and 


diffusion. These stages of gas breakdown are discussed in great detail in 


references 2-4 and 37-38. 


Ta Summary 


The cascade theory of laser induced gas breakdown has adequately 


explained most of the experimental results obtained to date, especially in 
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the case of high pressure and long laser pulses. The question of the 
effect of excitation on the breakdown process is not yet clear and 


further work needs to be done. 
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CHAPTER III 


DEFINITION OF GAS BREAKDOWN THRESHOLDS 


3.1 Outline 

Breakdown threshold has not been defined in the same manner by 
different workers in this field. These varying definitions as reported 
in the literature are discussed. Regardless of the definition used, 
however, the probability of breakdown is found to be a strong function 
of the laser power. In this work, the breakdown threshold, $¢ is de- 
fined as a 50% probability of visually observing the formation of the 
spark. This is the most common definition of breakdown threshold and 
by far the most convenient to apply experimentally. The effects of 
various experimental conditions on measured thresholds are considered 


and the procedures implemented to eliminate these problems are discussed. 


3.2 The Breakdown Threshold 

Breakdown is usually achieved by focusing the laser by a lens or 
mirror of short focal length to create a very high radiation field in 
the focal spot. The dimensions of the focal spot can be measured ex- 
perimentally or computed theoretically. However, there is sufficient 
inaccuracy in the spot size thus arrived at such that in breakdown work, 
agreement between two sets of measurements by different people within 
a factor of two may be considered very good agreement. This will account 


for the differences in the thresholds published. 
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The breakdown threshold, 4 is determined by the pulse length. 
For long pulses (> 10 ns) breakdown occurs when the rate of energy absorp- 
tion overcomes electron losses due to diffusion and recombination and 
energy losses such as elastic and inelastic coliteione. oe In this case 
¢ is usually given by the peak laser intensity in Fite For sufficiently 
short pulses, the losses are negligible and breakdown occurs when the pulse 
has sufficient energy to rapidly ionize the gas. Here, $ is measured 
in terms of energy flux ae is given in arene 

Experimentally, ¢ can be defined in several ways which may not be 
equivalent to each other. For example, breakdown threshold, 9 has been 
defined as the laser energy or power that will give a 50% probability of 
visually observing a eae oe Others have chosen other probability values 
such as 202°. Still others monitor the amount of light emitted by the 
plasma with photomultipliers. Breakdown can then be defined as the instant 
when a certain amount of light is emitted by the wikigneee Still another 
way of defining the occurrence of breakdown is the point when the total 
transmitted laser energy is a predetermined fraction of the incident laser 
energy’’®, Charge collection plates may be used to measure the total 
number of electron-ion pairs created in the focus with breakdown being 
defined as the certain value of the electron concentration. 

Regardless of the definition of breakdown threshold selected, the 
threshold measured experimentally is well defined. For example, Tomlinson 


and Damon measured the transmitted energy and plotted this against the 


incident laser energy. Their results are reproduced in Figure 3-L. ‘The 


sharpness of the curve near threshold indicates the sharp distinction 


between no breakdown and breakdown. 
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In the present work, breakdown threshold is defined as the laser 
power that will give a 50% probability of observing the appearance of the 
spark. $ was measured without preionization from very low \(25eTor::) 
to high pressure (600 Torr). The results are shown in Figure 3-2. Note 
that the threshold is a very strong function of the incident laser power, 
i.e., a small change in incident power will lead to a large change in 
the probability of observing a spark. In these measurements, a CO, laser 
was focused by a 5 cm focal length NaC& lens in Ar. Hacker, et. Hee also 
observed a change in probability from 0% to 100% with a 5% change in laser 
energy. In the present work, the laser is fired into the breakdown 
chamber and the number of sparks in every ten shots times ten is the 
probability of breakdown in percent. Because of the sensitivity of the 
probability of breakdown to incident laser power, in most experiments it 
is sufficient to determine the probability of observing a spark at two 
values of the laser energy and interpolating or extrapolating to 50% 


probability. Of course, these two probability values should not equal 02 


or 100%. 


3.3 Dependence of ¢ on Experimental Conditions 


Small amounts of impurities can alter the breakdown teehee 
These impurities may have low ionization potentials and can easily be 
ionized by the laser providing large amounts of initial electrons to start 
the avalanche. Avalanche breakdown assumes that there is at least one 


electron in the focal volume at the start of the laser pulse. Since the 
5 


~ 


=5.3 : 
focal volume is usually very small (V = 10 “cm’) while the background 


electron density n_ = ae ce without deliberate preionization, the 
fe) 
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probability of finding an electron at the focal volume to initiate the 
avalanche is small. Changing n, can change $ in a very complicated 
manner. Large particulates present in the focal volume may be ionized 
leading to formation of shock waves. In the region of the shock, 94 will 
be changed due to increased neutral atom density. 

These considerations indicate the need to condition the gas before 
measurements are made in order to ensure that the thresholds obtained are 
true avalanche thresholds and free from the effects of impurities, 
particulates and the statistical nature of breakdown at low initial elect- 
ron densities. Conditioning is achieved by firing the laser above threshold 
into the gas thus producing a spark consistently. The laser has to be 
fired into the breakdown chamber at a sufficiently rapid rate so that there 
is at least one electron at the focal volume at the start of the next pulse. 
The energy of the incident beam is then gradually attenuated until the 
desired probability of breakdown is reached. In gases conditioned in this 
way, the thresholds will be much more reproducible indicating the presence 
of priming electrons. That the conditioning is not due to dust or impurities 
that may have been blown off the lens or the walls of the chamber can 
easily be experimentally Verified. Breakdown measurements are usually 
performed in chambers with static gas fills. It is not unusual to observe 
that it is harder to break down fresh gas compared to gas that has been 
broken down several times. It is tempting to say that this is due to dust 
blown off from the walls of the chamber. If this is indeed the case, then 
we expect a continuous lowering of the discharge as we fire more and more 
shots into the gas. However, this is not observed. In actual fact, the 


threshold levels off after a few shots. Furthermore, if the laser is now 
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pulsed at a slow rate, the statistical nature of the breakdown returns 
and the gas becomes more difficult to break down. This suggests that 
conditioning is not due to dust or particulates blown from the chamber 
walls. A rate of 1 pulse/sec was found to be sufficient to condition 
the gas and assure the presence of priming electrons in the focus at the 
start of each pulse. Thus the thresholds obtained are true avalanche 


thresholds. 


3.4 Breakdown Threshold in Numerical Calculations 

In most theoretical work, breakdown is defined as the occurrence 
of a fixed electron density, Dd in the focus. The exact concentra- 
tion at breakdown is difficult to determine, however. Many workers 
just assume full ionization as a reasonable value. Charge collection 
experiments,on the other hand, indicate that 10-2 electron-ion pairs 
are created when a spark is produced. However, this may be different 
from the electron number at the focus at breakdown since it includes 
all charges created after breakdown. At best it sets only an upper 
limit to Ma? A charge collector was constructed in an attempt to mea- 
sure the electron density, n, as a function of the incident laser power. 
This will be described in more detail later on. However, the strong 
dependence of the breakdown probability on the laser power doomed 
this attempt from the start. Nevertheless, it managed to set a lower 
limit to Ma’ How this problem can be circumvented in numerical cal- 


culations will be discussed in Chapter VII. 
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Figure 3-1 Dependence of Transmitted Laser Energy 


on the Incident Laser Energy (7) 
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Figure 3-2 Dependence of Laser Energy on the Probability 


of Breakdown. Threshold is defined as 50% probability of 
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CHAPTER IV 


co, LASER SYSTEM AND ELECTRIC DISCHARGE PREIONIZER 


4.1 Outline 

The experimental setup used is described in detail. Parameters 
of the TEA co. laser system used in the breakdown measurements are given. 
Pulse shape, energy, mode pattern and stability of the laser output are 
also discussed. 

The electric discharge preionizer used in the experiment is studied 
in detail. Current and voltage across the discharge as a function of time 


are presented. 


4.2 Description of Experimental Setup 


The experimental setup consisted of three main components: a CO, 
oscillator-amplifier system, the electric discharge preionization source, 
and the breakdown chamber. A schematic diagram of the setup is shown in 
Figure 4-1. The output of the co, oscillator-amplifier system was focused 
by a 5 cm focal length NaC% lens into the breakdown chamber which had been 
filled with Ar gas at pressures varying from 100-6000 Torr. The laser was 
focused inside and outside the discharge region. The volume occupied by 
the discharge was ionized by the electric current and the gas outside it 
preionized by uv photons from the discharge. By varying the time delay 
between firing of the discharge and the laser and the distance between the 
focused beam and the discharge, a large range of initial electron densities 
Most of the measurements were made at 200, 400, and 600 


was obtained. 
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Torr. A small fraction of the incident beam was sampled by using an NaC 
beam splitter in conjunction with a photon drag detector and an energy 
detector which had been calibrated against a commercially available Gen-Tec 
energy detector. An attenuation cell placed between the laser and the 
breakdown chamber permitted regulation of the laser intensity incident on 
the breakdown chamber. Regulation was achieved by varying the pressure 


of propylene gas inside the eeu 


4.3 Pulsed Laser-Amplifier System 


The pulsed laser-amplifier system used has been discussed before in 


detail>. The laser oscillator was a helical TEA CO, discharge laser with 


2 
a peak power of approximately 1 MW and a total energy of 0.5 J. The 
output pulse had a rise time of approximately 60 ns and a width of 200 ns 
followed by a long tail several microseconds long. 

The amplifier was a Rogowski profile uv preionized TEA discharge 
unit. It had a typical peak small signal gain of 0.04 enue The oscillator 
output was made to pass three times through the amplifier providing a peak 
power gain of ~ 10 and an energy gain of ¥ 6. The initial part of the 
oscillator pulse was amplified more than the later part thus producing an 
output pulse with the same rise time as the oscillator pulse but with a 
half-width of only 100 ns. Oscillograms of the oscillator and amplifier 
outputs are shown in Figure 4-2 (a) and (b), respectively. Oscillator self- 
mode locking could not be observed in these oscillograms due to oscilloscope 
bandwidth limitation. 

Synchronization of the oscillator and amplifier discharge was 


important not only for maximum laser output but more importantly, to prevent 
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arcing in the amplifier discharge.” If the oscillator was fired before 
the amplifier, sufficient oscillator power was generated which can 
cause localized preionization and arcing in the amplifier. Synchroni- 
zation was achieved using multivibrator delay units. 

The laser output consisted of 3 to 4 transitions in the P branch: 
P16, P18, P20, and P22. However, most of the energy was concentrated 
in the P20 transition. The contribution of each line in the total out- 
put was not determined since it is not critical in gas breakdown 
measurements. 

The oscillator output energy varied less than 10 percent from 
shot to shot when fired at the rate of one pulse a second. The oscilla- 
tor-amplifier system was not as stable, the output varying by about 15% 
which is still an acceptable value. Figures 4-3 (a) and (b) show the 
variation in the output energy of the oscillator and the oscillator- 
amplifier system, respectively. 

Although the mode pattern of the output did not remain the same 
for all shots, most were "Gaussian-like" as observed on burn patterns 
on polaroid film. This was possible only with an iris inside the laser 
cavity which suppressed all other modes except the fundamental. A 
pinhole (< .5 mm diameter) and an energy detector was used to scan the 
beam energy distribution. Profiles obtained are presented in Figures 
4-4 (a) and (b). The half-width of the beam was approximately 1 cm and 
was not affected by the amplifier. 

By scanning the beam at two positions 1 meter apart, the beam 
divergence can be determined. The beam divergence was approximately 2 


mrad which is in close agreement with the theoretically computed value 


using Gaussian beam optics for the far field pattern. 
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4.4 Electric Discharge Preionizer 


The circuit diagram for the electric discharge preionizer is 
shown in Figure 4-5. A 0.2 uF capacitor charged to 25 kV was fired 
through a spark gap between two pointed tungsten electrodes 14 mm apart. 
The current in this circuit as measured with a current transformer was 
a damped ringing current with a period of approximately 200 ns. The 
current did not damp out sufficiently fast enough which interfered with 
the probe measurements to be described later. Hence, a 72 series re- 
sistor was added to critically damp the discharge current. With this 
modification the discharge current was measured to be less than 40 mA 
fifty microseconds after the discharge was fired. Peak current was 
2.8 kA and it was not sensitive to gas pressure. Oscillograms of the 
current are shown in Figure 4-6(a). 

The voltage across the discharge was measured using a divider 
consisting of a liquid resistor in series with a 1 2 carbon resistor with 
a divider ratio of 3,200:1. The divider was calibrated by applying 
the output of an EG & G pulse across it and the voltage output observed 
with an oscilloscope. The EG & G voltage was monitored using a 
Tektronix high voltage probe. The Tektronix high voltage probe could 
not be used to measure the discharge voltage directly because of induc- 
tance effects. Oscillograms of the discharge voltage are shown in Fig- 
ure 4-6(b). By integrating the product of the voltage and the current, 
the total energy dissipated in the discharge was estimated. 

The rapid release of energy in the discharge will cause an expand- 
ing shock wave to propagate through the gas. The shock will affect the 
breakdown ty duholds as they pass through the laser focus. It will be 


necessary to know the shock energy in order to relate the density in 


; ee pee «ie sae og ef i* ay #3 j<°) Dies bd BRgE Se | si nthe 3 > FY fit. at bo 


Ee Weetly aes we casa" ws 2 Atee thsi nna atic 


Hid favs brea fakhw dein ie + ebdaiak sia au Pon . ry 


~ 


eo eeteicg OC @. eoee .sapel perln wie i al sani 
aint. ha ie dns seznioekh ony quay, hi ty snag jy, ai ae 


AaeoQe gaia Qa l oa 62 Ss epanom aay iJ natdrhiy ‘a rei sn 8 é 
; . . ; iH nk “th 


‘ , ah, 


ad? 3o amare Piane .acyaaon) ayo aN seb ot #9 eit 4 
(aya * ong dy 


sah t » | f wiTe Le sé (Peat \ eR po vi ae _ my 


a9 
b- 
— 
és 
> 
— 
a 
= 
eo 

—- 


fy hw vee DRS L mitao ao dt o Ow ee eee “Senko have ee 


gaheey, ‘2 hothw és fen Bee eblih mii my sa yh Wy oka 


‘ 


& Sotau DSkoTkw AgwW SEI AOVou s bs 8 aif f seo a 


‘bide otoig spailie fold aiporsaeyl ach, ‘@ eer oxaabely ee ‘ id 
wahiad By sevkoud eiscevth apastoy adit pails ry aia o 


“Olt te OCs ote san} lou Measdeeth adn te 


, an Wis % tid ee any Love od4 Jue 34 atte 
(DRIES wae ewes iB: wl ah} 
~heeyes mm antes Tike oniadsekh wil) ok vee 
, 


40S Soatiy Lijit. ood afl (289 ED ayugsas 


/ sof aes SO ences. zanel ada ay wien nhs° 


the shock to the breakdown thresholds. 

The total energy converted to mechanical energy will not be equal 
to the total energy dissipated in the discharge but will be this energy minus 
the energy expended to ionize the gas and energy lost as radiation. This 
mechanical energy in turn can be related to the shock strength anywhere in 
the gas or to the gas density. If radiation losses are neglected then the 
total electrical energy converted to mechanical energy can be related to 
the total voltage across the discharge minus the voltage in the anode and 
cathode drops. 

The Pecumiqien employed to determine the total energy converted to 
shock energy is well known. Oscillograms of the voltage across the discharge 
were taken for the following electrode separations: 0.4, 1.4, and 2.4 cm 
at pressures of 200, 400, and 600 Torr. All other parameters were held 
constant (C = .2uF, L = luH, V = 25kV). At a fixed instant in time one 
could plot the dependence of the voltage on the gap size. These graphs 
are shown in Figures 4-7 (a) - (c). There is a linear dependence between 


the voltage and the gap size. 


Here, V is the voltage across the gap, Vi, (t) is the sum of the cathode 
g 
and anode fall voltages, E(t) is the average electric field in the plasma 


and & is the electrode separation. By extrapolating the data to =0 


we arrive at an estimate of the electrode voltage drops as a function of 


the total gap voltage. This relationship is shown in Figure 4-7 (d). Note 


that the cathode and the anode drops are independent of the gas pressure 


at the pressures considered. By subtracting this voltage from the total 
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gap voltage and integrating its product with the current, we arrive at 


values for the total electrical energy converted to shock wave energy. 
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The total energy dissipated in the cathode and anode fall regions was less 
than 1% of the total shock energy. On the other hand, approximately 85% 
of the initial energy stored in the capacitors was dissipated in the 7% 


series resistor. 
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Figure 4-3 Stability of Laser System Output 
(a) Oscillator Only. (b) Oscillator-Amplifier. 


The pulse rate is once a second. 
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Figure 4-4(a) Beam Profile of Laser Oscillator 
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Figure 4-4(b) Beam Profile of Laser Oscillator-Amplifier System 
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Figure 4-6(a) Electric Discharge Current.Vertical: 
1.8 kA/div; Horizontai: 200 ns/div (a) 200 Torr 


(b) 400 Torr (c) 600 Torr 
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Figure 4-6 (b) Electric Discharge Voltage.Vertical: 


25 kV/div; Horizontal: 200 ns/div (a) 200 Torr 


(b) 400 Torr (c) 600 Torr 
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Figure 4-7(a) Dependence of Discharge Voltage on Electrode Separation 
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Figure 4-7(b) Dependence of Discharge Voltage on Electrode Separation 
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CHAPTER V 
ELECTRON, ION, AND NEUTRAL ATOM DENSITY: 


THE REGION OUTSIDE THE DISCHARGE 


5.1 Outline 

Diagnostics performed to measure the electron-ion density in the 
region outside the electric discharge were performed. Probe techniques 
were used as these afforded good space and time resolution. 

A probe model is developed which extends the theories reported in 
the literature. The theory applies to high pressure plasmas flowing 
to a cylindrical probe surrounded by a thick sheath. 

The strength of the shock wave created by the discharge was measured 
using a piezoelectric pressure probe. Using shock theory, the neutral atom 


density at the shock can be calculated. 


5.2 Introduction 


Recent measurements of uv photoionization-created plasmas have 


oe 
employed microwave interferometry or charge collection >”, These methods 


do not give good space or time resolution which is important if we want 
to accurately determine the initial density in the focal volume that will 
initiate avalanche laser breakdown. These measurements indicate that uv 


i Avs i h < 10" me Electron- 
photoplasmas typically have densities in the range n £ cm. Cc 


ion densities as low as 10° ge can be measured using swinging Langmuir 


3 
probes with good space and time resolution . 
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5.3 Probe Theory 


Probe theories have been developed for a large range of ion densities 
in high pressure flowing pivenae o". A review of these theories is avail- 
able in the literature?. in high pressure probe theories, parameters 
such as probe geometry, plasma flow, bias voltage, sheath thickness, 


hydrodynamic and thermal effects, and so on should be considered. These 


parameters are discussed in terms of the following nondimensional quantities: 
R. = v_L/ (ukT/e) = electric Reynold's number 


and yy = i ys 


where 
vee plasma flow velocity 
L = characteristic length = ai for cylindrical probes 
Ms = probe radius 
u = ion mobility 
K = Boltzmann constant 
T, = electron temperature 
e = electronic charge 
dy = (KI fue) = Debye length 
ar permittivity of free space 
n = ion density 
es = probe bias voltage (negative for ion collection) 


For low densities such as in uv photoplasmas, the current to a 


moving probe is described by the so-called sheath/convection theory which 
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is applicable when Roa 7>L. In this theory it is assumed that a thick 
sheath is established around the probe (sheath thickness cee 2 The 
sheath is thick if Re aoe ees ee 

The current to a cylindrical probe may be limited by the space 
charge or by the electric field in the plasma sheath. In the first case, 
R a/ Vx<< 1 and the sheath electric fields are so large that ion motion 
in the sheath is radial and all ions that enter the sheath eventually reach 
the probe. The case has been studied by Clements and phy In the second 
case, Ral Vx>> 1 and only ions that enter the sheath within a limited range 
of impact parameters can reach the probe, the rest passing through the 
sheath. This is shown in Figure 5-1. In the transition region, Ro/ Vx ae 
it has been suggested that one should use the model predicting the smaller 
probe current or equivalently, for a measured probe current, the theory 
predicting the smaller ion density may be ubeeee 

The space charge limited and field limited models are very similar 
and a closer look reveals that the two are not inherently different but are 
limits at high and low densities. The following is an attempt to extend 


these theories and account for the ion current to the probe for all values 


of Ra/Vx. 


5.3.1 Derivation of Formulae 


3} 
The same assumptions made by Smy and Noor are made, namely: 


1 The current is convection dominated, i.e., the current due 


to diffusion of ions into the sheath is negligible (Roo Yared). 
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Zs The sheath is thick and hence, hydrodynamic effects are 
ete =I 22 
negl b 
gligible (ro >> ri Ro Qey see Vu 
a The sheath is cylindrical with uniform charge density en and 


there are no end effects. This assumption is most reasonable 
when R,o/vy >>1, i.e., when n is very small. When Ralvyx << us 
ye has shown that the sheath charge density is not uniform 
but increases from a small value near the probe to a maximum 
before decreasing to zero at the sheath edge, r = to" There- 
fore, the greatest error in the model should be expected when 
n is large. 

4, The electric field and potential at the sheath edge are zero. 
In the case of large bias voltages which are usually the case 


in experiments, the quasi-neutral transitional zone where there 


is a voltage ~ KT /e°?’ can be neglected. 
The charge per unit length, -Q on the probe is: 
2 
Or me(r.” i te, (Sh) 


where r is the radius at the sheath edge and os is the probe radius. 


The electric field per unit length due to this line charge and the surround- 


ing plasma sheath is: 


Ee=- ave 2 He TR AS 8 (5.2) 


Integrating, the probe voltage becomes: 
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2 Be 
ner ty Z to l i* 
V —— SS Se — 
pie iomtenory (yu oe @t sey a (5.3) 
p Pp ig 
re) 
This equation can easily be solved for ros by iteration. Figure 5-2 


shows the dependence of r /r_ on 2e V /ner 
OF Dp op 


> 
From the electric field E and the flow velocity v the ion orbits 


iE 


inside the sheath can be determined. The velocity inside the sheath is: 


> > ie 
v=v.+uE (Sea) 


while outside v = Ve: Taking the x-axis to be parallel to Veo then 
in cylindrical coordinates: 

> > _. + PATI 

Ver= Vets v,_ (1,.cosé 1p sin6). (575) 
> 
E is radial and hence, 

2 2 
ei ne re) 5D TF 
E=- oem aaa eae ) 1: (5.6) 
fe) 
Therefore, 
2 2 
as eS neu foe é } 
Megs i _[v,cos® - 2a Ser - iy v sind CHIP) 


This equation can be converted to a non-dimensional equation by the trans- 
ner 


formation r =p) (geek Separating the r and 6 components and 
oof 


converting to non-dimensional variables: 
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2 2 
ao ( foe ) [v,cose - ( ct: e ( Poe )] ( 

- a Vv st 5.8a 
dt ner * u £ 2€ V¢ iB 0 ) 
de _ 2€ Ve v,-siné 
dt 2 0 6 38b) 

ner u 


Dividing the first equation by the second yields the differential equation 


for the ion orbits: 


2 


dp pe 3 
46 ‘Tp ecOto, + ‘ 5 csc@ = cscé@ (5.9) 
fe) 


where P, is the dimensionless radius at the sheath edge. 


3.3.2 Asymptotic Solutions 


The field limited case occurs when n> 0, Eerek he tond \eote 
This is the case of very weak shielding by the plasma sheath. When 


pe © the third term on the left of Equation 5.9 approaches zero. Its 
) 


solution is: 


osinéd = 6 + constant (516) 
and 
2 Fo Y 
1/2 = © (ner_“)(— ) (5.11) 
oe P TD 


3 
This is the field limited solution of Smy and Noor . 
The space charge limited case (sheath/convection) is described by 


69 +0, n>, r +0. The asymptotic solution can easily be obtained if 
Oo 
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we solve for d0/dp: 
zu 
Ce pa siné 
on Se eee ee ec a, ee, C5125 


é 2 
P, Pcot® + cseb(p + p62) 


Taking the limit as pass 0, 


CR AO eu 10 (5.13) 


constant (5.14) 


@ 
i} 


I/2 = 2nev Ba ( et ) (5.15) 


This is approximately the same result as in reference 6. 


Deo.3 Transition Region 


The full differential equation (Equation 5.9) has to be solved in 
the transition region. This has been done numerically using a fourth order 
Runge-Kutta method. There are two classes of ion orbits. Depending on the 


initial value of the impact parameter Y = p sine , the ions either reach 


the probe or pass on to X = pcos® = ©. The solution which separates these 


two classes of orbits enters the sheath at the point (05°81) where 
P51 70° It is easy to show that this orbit should include the point 
(p ir? which is the solution of the equations r=0 and 6=0. The radial 
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motion of the ions is due to the net radial velocity under the influence 
of the electric field and the fluid flow. The ions move towards the probe 
as long as r is negative. If ft becomes positive for 0 < 6 < 7, then 
the ion will never reach the probe. The maximum radial component due to 
the fluid flow is equal to the flow velocity itself and this occurs when 
6=0. Clearly, the condition r=0, 6=0 gives the ion orbit which separates 
the orbits of captured and uncaptured ions. 

In the actual computer program, the solution was started at Peeie 
and 0=0 in where ha was a very small angle (= nore radian) which will 


not cause computer overflow. The solution was then generated by going back- 


wards in time. This yielded (0912 %)> The values of 0 and 


01 
; ef : F = 
Y917P 91 Sin® 91 =P 98in8 54 are plotted in Figure 5-3 as a function of Po - 
For small values of ee You + a1 which is the field limited solution. 


When i is large, You * Po and 81 + 1/2 which is the space charge limited 


case. From the values of Y in Figure 5-3, the ion current to the probe 


ol 


can now be written: 


I uner 
Q g Yo1 ( 2EnV ) Bees 
yY L Z 
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fe) Pp 


Thus, for particular values of hey yw and ae Fquation (5..3) can be used 


to solve for r /r. as a function of n or eas can be read off from 
es Leny. 
; One ; Rs 
Figure 5-2. Knowing 9, = ner ou and You (see Figure 5-3), plots of I/% 


versus n can be constructed. From these plots n can be read off for 
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each experimental value of I/2. An example of such a plot is shown in 


Figure 5-4 together with I/2 predicted by the space 
field limited theories. It is to be noted that the 


to all values of R,a/Vy. In the transition region, 


charge limited and 
present model applies 


the difference between 


the predictions of the present model and previous theories mentioned above 


is approximately 50%. Figures 5-5(a) - (b) show probe characteristics for 


two values of the electron density. The agreement with experiment is 


quite good. 


The sensitivity of the probe currents calculated numerically to the 


exact form of the sheath density distribution has been tested by trying 


other charge density profiles. Specifically, the following profiles were 


used. 
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Pcharge ae _ % 
> 1/2 iE ok F=1. .5 
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Coe Le) 
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me (Se17£) 


These profiles cover a broad range of possible charge distributions in 


the sheath including cases where there is charge depletion and no charge 


fepietions. 


In all cases considered there was no significant change in 


the predicted probe current. As expected the largest error (~ 20%) occurred 


at the high density limit. The error may be due to the fact that in this 


limit r /r < 2 and hence, the condition of a very thick sheath is not 
op 
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Satisfied. A possible solution to this problem is changing the assumed 
constant flow to the following inviscid solution which takes into account 
the effect of the finite size of the probe. 


os 


ES <2 Zee : 
eee habs dag h /r-)cosé - tte + r*/r)sine } (5515) 


This was not done in this work since the measured probe currents were 


far from this asymptotic limit. 


5.4 UV Photoion Measurements 
5.4.1 Experimental Setup 

The experimental setup used is depicted in Figure 5-6. A serious 
problem with using probes to measure uv photoionization is photoelectric 
emission from the probe if it happens to be exposed to the same uv radiation. 
To overcome this problem, a metal plate was placed between the discharge 
and the probe. The plate blocked all radiation from the discharge except 
that which passed through a 3.2 mm diameter hole. The rotating probe was 
kept in the shadow of the plateup to 50us after firing the discharges. At 
this time, the discharge current was less than 40 mA. This low current 
minimized the possibility of unwanted photoelectric emission from the probe. 
A volume ionization in the shape of a cylinder was produced in the measure- 


ment chamber. The photoion density was measured 1 cm from the metal plate. 


5.4.2 Probe Design 


The probe had to be sturdily constructed and well balanced to withstand 


the tremendous stresses on it when operating at very high rotational speeds. 


The probe was driven by an electric motor outside the measurement chamber 
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at about 7200 rpm giving a probe speed of approximately 100 m/s. Oil 
lubricated rubber seals ensured that the measurement chamber was vacuum 
tight. The probe can be operated safely up to 10,000 rpm. At higher 
speeds, the oil lubricated seals started to fail due to the excessive 

heat from friction. The probe was made of 0.025 cm diameter stainless 
steel wire. This was placed and cemented inside a ceramic-lined stainless 
Steel tube kept at ground potential. The tube also served as an electro- 
magnetic shield against noise. Figure 5-7 shows the probe circuit and the 
probe construction. A 55 volt battery provided the negative bias and the 
total resistive load was varied from 6002 to 10k2. Electrical connection 
between moving and fixed parts of the probe circuit was through a mercury 
pool. Measured probe currents ranged from approximately 250nA to 20uA. 

The time resolution of a moving probe is approximately a few times 
the time it takes the probe to traverse a distance equal to its diameter 
For ty = 1.27 x om cm and v = 10° cm/s, this is about 2.54 x ene ey 
Hence, with a probe time resolution of = 5yus, the RC time of = lus for the 
probe circuit is adequate. 

Figure 5-8 shows typical oscillograms of the probe current at 
different times after firing the discharge. The upper traces are probe 
current signals while the lower traces are reference signals obtained from 
a phototransistor when a He-Ne beam was chopped by the probe arm. This 


reference signal located the probe relative to the plasma. 


5.4.3 UV Photoion Density in Ar 


Since the probe length was greater than the diameter of the cylindrical 


plasma, the effective length of the probe, Lee (see Equation 5.16) at any 


time when the probe is in the plasma is equal to a chord of a circle as 
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shown in Figure 5-9(a). For convenience, this length was estimated by the 


distance, 2, corresponding to the full width at half maximum of the probe 


Signal trace. Hence, 
Q = 2 = vAt C5: 49)) 


where At=t,-t, is the FWHM of the signal trace (see Figure 5.9(b)) and 

v is the probe velocity. Referring to Figure 5-9(c), let Ry and R, be the 
radii of the plasma at times t, and t, such that At=t,-t,- At ty and ty 
the probe is positioned a distance Ty and Tos respectively from the center 
of the plasma as shown in the figure. Then 2 = r,tr,: Since TS f oe Sho 
and since the plasma does not expand rapidly at the time the measurements 


were taken (t>50us) (see Figure 5-1), then it can be seen that 2% reason- 


ably estimates £L gge this value of & will yield the ion density, Daye (t) 


averaged over a diameter. 


Probe currents as a function of position and time can also be obtained 
from the oscillograms. From these data, the density as a function of 
position and time can be obtained after an Abel inversion. However, the 
procedure described above was adequate because the results were not suf- 


ficiently accurate to warrant an Abel inversion. 
The plasma expanded radially due to ambipolar diffusion, & in- 


creasing in time as e (see Figure 5-10) which is to be expected for an 


unbounded gas diffusing in two dimensions. 
It will be shown later that T, *’ 1 eV and using the experimental 
-2 2 
= 55 volts, v,+10"cm/s, ste cig cm, up =2 cm /V-s and 


1 aa 3 


parameters V 
10 - 
a x <1 yield ns5xl0 cms. 


Dae = 
p=100 Torr, the conditions R a x >1 and Ro 
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Therefore, the theory discussed in Sec. 5.3 applies when n < yO: ema: 

Measurements were made for pressures of 100, 200, 400, and 600 Torr 
at distances of 2.5, 3.8, 5.1, 7.6 and 10.2cm from the discharge. The 
discharge chamber was not isolated from the measurement chamber and, there- 
fore, the pressure in these two chambers change together. Thus, it is not 
expected that the source spectrum will be the same for the different 
pressures considered. Impurities such as mercury from the probe connector 
and oil from the seals could also leak into the measurement chamber. No 
attempt was made to eliminate these impurities nor were attempts made to 
determine to what extent these impurities affected the results. The reason 
is because the breakdown measurements reported here were done in systems 
which were far from pure. 

The results of the measurements are shown in Figures 5-11 (a) - (d). 
The points are experimental values while the solid curve is a fitting of a 


rans 5 
theoretical expression to be discussed later. The ion mobility was 


chosen to be: 


Loe 20 7 AI (5220) 
The densities obtained are of the same order of magnitude as those obtained 
by Seguin, et. Hee in a laser mix. The results of Babcock, et. wie are 
about three orders of magnitude smaller. The differences could be due to 
different current pulse magnitude and shape, total discharge energy, dis- 
charge gases, distance from source, target gas, and the gas through which 


All these parameters affect the emitted uv spectrum 
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the uv photons pass. 


and intensity and the photoionization efficiency 
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5.4.4. Curve Fitting and Extrapolation 


The photoplasma decays due to recombination and ambipolar diffusion. 
Diffusion along the length of the cylindrical plasma can be neglected 
compared to radial diffusion or recombination losses. The ion density 


decay is described by the following equation: 


(t) (5221) 


where A. is the radial diffusion length, Dy the ambipolar diffusion 
coefficient, and Op is the recombination coefficient. Equation (5.21) 
can be solved if an expression approximating We can be found. 

To do this consider a two dimensional, freely expanding plasma, 


where diffusion is the only loss mechanism. Then the density distribution 


is eeeegant” 
n(r,t) = —7-“—  exp{- rey (5.22) 


vee: ’ : 
R is the plasma radius at t=0 and R =R +4D,t is the mean square radius, 
fe) 
i.e., the radius where n equals 1/e of its central value. Hence, we 


identify p7eR?. From Figure 5-10 which gives 2% as a function of time, 


Drointeands 2 dhe Ss 
D, can be calculated using oe = R a +D,t. Ro = .22cm is just the 


A 
radius of the uv port. Assuming that 1,=300K"”, the electron temp- 


erature, T can be estimated using 
e 


kT. ‘i 
. i es (5.23) 
Dy ui, ¢ = ya Gal me ) 


The values of T consistent with the data in Figures 5-11(a) - (d) can 
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be found in the inserts on the figures. 


Averaging n(r,t) over a diameter: 


D 2/2 f R 7 co 
Nave (t) ey f n(r,t)dr ~ os fn(rst)dr = R ie n(r>t)dr 
fo) 
(See) 
or 
n (0) 
_ ave 
Nive t) 7 D, (5425) 
wae 
R 
fe) 
By assumption, this satisfies: 
“Pave od oie 2 ped Pallave St (5.26) 
dt A 2 ave 1 oR 2 AD (t)) : 
Ne 4 fo) A 
Therefore, 
2 i 2 Lae 
N ie (R, + 4D,t) vag Q C527) 


This is the desired expression for the radial diffusion length. The 


solution of Equation (5.21) using Equation (5.27) is: 


SS ——————— S525 
B agent? ( ) 


4D,t R ADrt 


A fe) 
(1 + ——;) [a ,n (0) —— &n(1 + ) +1] 
R Z R ave 4D, R 2 


oO 10) 


By a suitable choice of Op and Nave (9)> Equation (5.28) can be 


fitted to the experimental points. The solid lines in Figures 5-1l (a) - 


(d) were obtained using such a procedure. There is reasonable agreement 
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between experiment and theory. The values of Op used varied from 
-8 =f nS 
1.0x10 to 2.7x10 cm/s. These are typical recombination coefficients 


for dissociative molecular recombination prOceosS” 


Peco error in 71 a (t) 


The uncertainty in 9) will be due to errors in oe 
error in Has error in estimating 2, and variations in the current from 
shot to shot ( 10%). There is error in estimating 2 because it is not 
exactly equal to hee as discussed previously. Taking a relative error 
of 10% for u, 20% for 2 and taking the field limited expression for the 
density for example, the uncertainty in Nove can be estimated. In the 


-.//7 -.46 An 
¢ 


3 A A 
field limited case n ) ie . Hence, cane a7 “+ Ss, 


Ar 


+ .46 aa wa, a 35%. This is typical for probe measurements. 


5.4.6 UV Photoionization Mechanism 
Figure 5-12 shows the initial photoion density obtained by extra- 
polation of the data to t=0 using Equation (5.28), plotted against the 
distance from the discharge. Two photoionization mechanisms can be con- 
sidered: single photon ionization of Ar or of impurities and two photon 
ionization. For single photon ionization of Ar or of impurities, n 


a -upd : 
should vary with distance, d, as “5 e@ UPC where p is the pressure of 


d 
the uv absorber and yp is the absorption coefficient. No single value of 


ave 


u will fit the data of Figure 5-12 for all pressures considered and, 
therefore, we conclude that Ar is not photoionizaed by a single photon. 
On the other hand, impurities cannot have the large pressure required to 


fit the data. Hence, single photon ionization is ruled out. 


ee - Le * ie | 
; A st Bi oer ‘ - a 
Sy ae ee i PDL 
ak ‘ 4 a : Va; 
ry i * 
ONT 
nf 4 


(ROL) 7 at o7ors of Suh ad ‘kkhe 4, waa ae sai 
mee ame Ttied salts. 22 anosasbtay bae ‘ea snksnmbie: - pes 
Out at pr seus oat 2 goaitenites fk Bom et seat hos ) 
soux8 avitelet » anitdalr  .¢herogvaiq belamuys sae, ay “iat co 

‘ead tot eoledsiges hettotd bie ad okt yoked fas , 3007 ’ 


Sty Al «ieee 89 BBD w at eda kas gen ‘ofa aig 


E peat , 
. * af 
| ' rh. FY ; 4:\ 7, Mi me on: 
ee! eee ae = , orev fom ' ay, ae 
| " | 


esoaeoiuessm odew] int Iackews at ghar. We, 


gesaednot 

Ny “oartgete A benistdo voiengh meko: key fetann yet ‘wade * 
, nate yemkege heagolq , (48.2) notseupe gotow Const om stab 

| ers ot ve ates bast es moktemtrogysotty owt sqonidaath at, 
| nese ows bits eoltiung mr tu yo TE to sottaxkuos aes 


A old dzoget ts to “A 36 nok 4tes bow ct soso ny of ghka ‘sot 


a 


bi % eriiwe9>9 ott eh q every “i, - ei ab erateteit ) hat 
way ’ B, 
ule Be ater ba dial OW, »IeekolTieds eet aia ad EY at: ete > 
a i L oat betaBbeno> anmueas: ty Tle 104, SIAR onvabt to: 
¥ : 
i ek tans A ed — noLos od Jor wh oh yoda 
Pe ay BN 


| os atime 26 auyel sd Mier 2 crue) eukshy =i . 


ion a dxrx ef od qu Each. ose,» 
Niet 


oo ‘, : 


For a two photon process, the rate ions are created, W, is proport- 


ional to the square of the bight flix, 


2 
W = CF (5129 ) 
and 
dn 
——\= Wn 
dt gas C530) 


From the conservation of energy, 


dF ZEe 
dr y nhs ie W Comoe) 


which is applicable for a point source of photons. The solution to this 


set of equations is: 


2 
-1 
F(r) = F [((1+2n. Cr F) => - 2n. CF r (—)] (r32) 
O gas oo 2 Sass e077 2 
2 eo ig -2 
ee Gr Coated Bor aan Oe 0), be - see See ( zr, )] (S33) 
6) 


where ris the source radius (spherical), F=f (r)) and At, is the dura- 
O 


tion of the discharge current. Since Ar is transparent to uv longer 


i 2 
than 800 A, 2n_ CF r <<l. Hence, 
gas oO 


ee Sond GS (5034) 
n =CPr r Grae td ur 


The solid lines in Figure 5-12 is a fitting of Equation (5.34) to the 
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-4 
data. It can be seen that an r dependence is not incompatible with 


the data. 


The effective photons may come from the Ar continuum around 1200 


23524 
fi - . When a CaF, sindow was placed on the uv port, the probe current 


was reduced to 10% of its value without windows. However, it is known 


that a high uv flux such as that near a discharge degrades the transmis-— 


sion characteristics of uv windows such as Gal aih It is then very 


probable that the actual transmission limit of the window used in the 
experiment was higher than that published for CaF, ( =1100 aye This 
could explain the reduced probe current with the window in the uv port. 
D lasers, it was at first thought that the 
effective photons were in the region A<1100 A since the arc-suppressing 
WEP As 


In fact, in uv preionized CO 


effect of the uv discharge was inhibited by CaF, windows 

The fact that 10% probe current still persisted even with a window 
precluded the possibility of the ions coming from the discharge direct- 
ly via an ionization wave. When a NaC window was used, no ionization 
was observed. This means that A<1700 A. 

All this discussion is only speculation, however, A more detailed 
study of photoionization in Ar has to be performed before the question 
can be definitely answered. In breakdown measurements, it suffices to 


know the ion density without regard to how they are created. 


5.5 Pressure Probe Measurement of Shock Wave Strength 


5.5.1 Pressure Probe Characteristics 


The shock wave created by the rapid introduction of energy into 


the discharge column was measured using an LD-25 blast wave pressure 
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transducer. The transducer had been calibrated in a shock ye and 


the response was within the manufacturer's specifications of .15V/psi 
ec 2.2 V/atm. 

Figure 5-13 shows a typical pressure transducer signal response 
to a weak pressure jump. The oscillations are due to vibration of the 
piezoelectric crystal at its natural frequency after it is excited by 
the shock. 

Measurements were made with the piezoelectric crystal (lead zir- 
conate titanate) face perpendicular to the wave direction. Hence, 
recorded pressure jumps were the sum of the actual pressure jump plus 
the equivalent pressure exerted by the fluid flow. | 

The linearity of the transducer was checked by plotting the trans- 
ducer signal against the initial energy stored in the discharge capaci- 
tor bank. Assuming the energy in the shock remains a fixed fraction 


of the initial charging energy, Figure 5-14 shows that the transducer 


had a linear response. 


5.5.2 Measurement of Shock Strength 


The shock wave was monitored for pressures of 200, 400, and 600 


Torr in Ar at distances varying from 0.8 cm to 10 cm from the discharge. 


The total transducer signal, AP, is: 


2 
A u 
ee AP 4 Fil (5.35) 
Py a Py 


where Ap is the actual pressure jump, Py is the ambient pressure, Py 
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is the mass density of the compressed gas and u is the gas velocity. 


For a perfect gas with constant specific heats the following relations 


hoid: =) 


Py (y4tL) e,- (y-1) Po 


ls fh tEAM ED Ne (5.36) 
0 (y#l)e- (y-1)P, 
Py (y#L) (P,/P9) + (YH) 
B,  C-1)@,/P)) + OFT) eee 
Meee Ruy = (Pepe yee) joe (5.38) 
OF. LO" p P4 


0 


where, P and PS? Po» u_ are the pressures, mass densities, and 


u 
Mee? 0 


velocity of the gases with respect to the discontinuity of the compressed 
and undisturbed gases, respectively. y is the ratio of the specific 


heats. From Equation (5.38): 


Py“ 


p p 
Sey (poe) ieee (5.39) 
Po Po oOterat oF, £0 


where Ap = P1-Po 


Hence, 


Poe glee: 


Fo “fi Po 


(5.40) 
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Solving Equations(5.40) and (5.41) yields: 


A 
ie. (5.42) 
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Equation (5.42) gives the actual pressure jump in the shock from pressure 
transducer signals. The results are shown in Figures 5-15 (a) - (c). 

Theoretical investigations have been performed for spherical shocks 
by Brode and by Pigoseee for line sources. For spherical shocks, 


Brode finds: 


Poe aS 7h W110 1269S oa ae 
aire Ae ‘ 
where i = r/e and oe = E/p 9: E is the total energy converted to shock 


energy. This has been determined previously. (Chapter IV). 


For cylindrical shocks, Plooster reports: 


por) 2 ofa (sent? ® - ay} (5.44) 
Y Pg 
where 
: eee o8/5, gue 8/5 é SEY Mae r/R, 
12 
Ry = (E/ YP 9) 
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Plooster suggests c=.7 and b=3.94 for air. The same values of c and 
b have been used here for Ar. 

The solid lines in Figures 5-15 (a) - (c) are from Equation (5.43) 
and the dashed lines are from Equation (5.44). It is to be noted that 
the shock seems to start out as a cylindrical shock and becomes spheri- 
cal at large distances. This is reasonable to expect since the discharge 
looks like a line source at close distances but looks like a point source 


at large distances. 


5.5.3 Discussion 

The sudden increase in gas pressure when the shock front arrives 
at the laser focus should reduce the laser-induced breakdown threshold. 
This should happen because the pressure jump due to shock corresponds 
to an increase in the gas density as can be seen from Equation (5.41). 
The dependence of the breakdown threshold on the neutral density is 
known or measureable (see Section (2.5)). Therefore, it may be pos- 
sible to measure rapidly changing neutral atom densities with a high 
degree of spatial resolution with lasers. This indeed is the case and 


this will be discussed in greater detail in later chapters. 
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Figure 5-1 Ion Orbits in Plasma Sheath 


(a) Sheath/Convection Regime 


(b) Field Limited Case 
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Figure 5-5(a) Probe Current Versus Probe Bias Voltage. 
Dashed line is from Equation (5.16). The solid line and 


the experimental points are from Reference 3. 
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Figure 5-5(b) Probe Current per Unit Length as a Function 


of Bias Voltage. Solid line is from Equation (5.16). 
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Figure 5-8 Photoplasma Probe Currents.Typical Oscillogram 
of probe currents at two different times. The lower 
traces are the reference signals (see text). Sweep 


speed is 20yus/div. 
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Figure 5-1l(a) Photoplasma Density Versus Time 
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Figure 5-11(b) Photoplasma Density Versus Time 
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Figure 5-1l(c) Photoplasma Density Versus Time 
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Figure 5-11(d) Photoplasma Density Versus Time 
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Figure 5-12 Variation of Photoplasma Density with Distance from UV Source 
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Figure 5-15(a) Shock Strength as a Function of Distance. 
The solid line is from Reference 28 (spherical shock wave) 


while the dashed line is from Reference 29 (cylindrical shock). 
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Figure 5-15(b) Shock Strength Versus Distance. The solid 


line is the theoretical prediction assuming spherical shock. 
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Figure 5-15(c) Shock Strength Versus Distance. The solid 


line is the theoretical prediction assuming a spherical shock. 
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CHAPTER VI 


ELECTRON, ION, AND NEUTRAL DENSITY: IN THE DISCHARGE 


6.1 Outline 

Interferometric measurements using a CO, laser (10.6um) and a 
He-Cd laser (441.6nm) are described. A Mach-Zehnder interferometer 
was used to measure the electron density at early times. A Jamin inter- 
ferometer was used in conjunction with a He-Cd laser to measure the 
neutral gas density. Fringe shifts in the latter case are contributed 
by electrons, ions, and neutrals. 

Absorption measurements were performed to estimate the discharge 
electron temperature. The absorption mechanism at 10.6um is inverse 
bremsstrahlung which is related directly to the electron temperature. 

Finally, a rotating Langmuir probe measurement of the electron 
density for late times is described. Interferometric and probe results 


are complementary. The plasma decay mechanism is discussed. 


6.2 Interferometer Principles 


An interferometer measures the index of refraction of a given 
transparent sample. All interferometers operate on the same basic prin- 
ciple. A light beam is divided into component beams (wavefront division, 
amplitude division, or polarization Aivislon tend then recombined. The 
sample is placed in one of the component beams. Interference fringes 
formed by the recombined beams are observed; changes in the optical re- 


fractivity of the sample are seen as shifts in the interference fringes. 
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6.3 Interferometric Measurement of Electron Densit 
_ ee eeemennte Ot EtLectron Density 
623. lL Optical Refractivity of an Ar Plasma 


The total index of refraction of a plasma is the sum of the indices 
of refraction due to electrons, ions, and heltralsas For an argon plasma, 


the index of refraction is given be 


p=l+ (At B/°) n - cx n (6.1) 

a 
where 

A = 1.039 x 10 “em 

Papas co x 10° em 

ne = 1 + 0.8 n, 

a n ion 

=n +0.8n 
n 

n = electron density 

ed neutral density 

Cc = 4.48 x lOesecn 

X} = wavelength of light source (cm) 


The contribution of the charged and the neutral particles in a 
plasma are strongly dependent on the wavelength. Hence, by using two 
light sources at different wavelengths it is possible to separate the 
contribution of the charged particles from that of the neutrals. 


be -3 
With ST = 4.416 x10 oer (He-Cd laser) and do = 1, 06 x 10> cm 


(CO, laser) the indices of refraction are respectively: 
- -23 
Hee betty 009 x LO ve he a7 SODEXO n (6.2) 
i n 
ly = 1 + 1.039 x 100s ae oy Shea hie te (6.3) 
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The contribution of the electrons to the total index of refraction at 
10.6um is three orders of magnitude greater than that of the neutrals. 
Hence, if the change in the neutral density, An? is not much greater 
than the change in the electron density, An, the change in the index of 


refraction, Aun» is: 


Wy, = 1 - 5.033 x 1077 in 


Z (6. 3a) 


6.3.2 Mach-Zehnder Interferometer 

Figure 6-1 shows the setup of a Mach-Zehnder interferometer. Coated 
NaC beam splitters were used for M, and M,. The discharge was placed in 
an arm of the interferometer. Scanning of the plasma was made possible 
by placing the whole discharge system on a movable base. Changes in the 
index of refraction of the plasma caused shifts in the fringe pattern 
which were observed with a Ge-Au detector. A filter provided rejection of 
the discharge emission. The number of fringe shifts P, is: 

Au,’ 


ee (6.4) 
az ho 


where 2 is the optical path length and ro is the wavelength. One 


fringe shift implies Oo = 10.6um) 


-2 
<pnios = 2) xX TO. Cen (Gon 


6.3.3 Cylindrical Geometry 


The plasma was assumed to be cylindrically symmetric and since it 


was illuminated "side on", the fringes observed as a function of the radius 
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had to be subjected to an Abel inversion to obtain the refractive index 
as a function of the radius. Referring to Equation (6.4) and Figure 


6-2, the fringe shift at r =r. is: 
0 


an s 
Z 0 
D re Au, (x) rdr 
rics ee oer (6055 
2 a vee pe 2 
i 


This is Abel's integral equation. The solution for this eduatilonw isi: 


ape Be ae semen eee 
Au, (r;) = met dz (r a ) dr (6.6) 


From Equation (6.6) and Equation (6.3a) the radial electron density 


distribution is obtained. 


6.3.4 Electron Density Distributions 


Figure 6-3 shows a typical oscillogram of the fringe shifts. The 
fringes are counted from the right, _To avoid confusion regarding the fringe 
turnover point, fringes were counted no earlier than 30us after the dis- 


charge was fired. Figures 6-4(a)-(c) show the fringe shifts plotted 


against the radius at four different times and pressures of 200, 400, 


600 Torr. A polynowpial of the form: 


P, = are o pre +c (6.7) 
was fitted to the data. This is just a fourth order power series expansion 


with the coefficients of the odd powers of r equated to zero since they 


lead to a discontinuity in the derivative of nat r-=0. This is clearly 
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physically impossible. Substituting Equation (6.7)into Fquation (6.6): 


Ly 2 2 


: vat x (a'r - b') (6.8) 


is 
Au, (r) = -6.7 x 10 (r. ) 
where a' and b' are expressed in terms of as Lbs eand Ts 
The results are depicted in Figures 6-4 (q) - (£). These are 


comparable to the results of Offenberger, et. ae ina co, laser spark 


on AY. 


6.4 Absorption of 10.6ym Radiation 


By blocking the reference arm of the interferometer the absorption 
of CW co, radiation by the discharge plasma was measured. The absorption 


process is inverse bremsstrahlung and the absorption coefficient, a, at 


F > , 
where n is the electron density in cm , Ty is the electron temperature 


in eV and « is the photon energy in eV. The transmission of the laser beam: 
af 


T = exp[-Sa(2)d2] 


allows one to compute the electron temperature since the density is known. 


Figure 6-5 is an oscillogram of the transmitted CW signal showing absorption. 


By using a large collecting lens (£-= 20 fem,cf/3.43)) the effect. of refraction 


at the higher electron densities was reduced. This was verified when the 


beam was apertured and the percentage transmission remained approximately 
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the same. 

Figure 6-6 shows the measured transmission at the pressures 
considered. The averaged central electron temperatures were obtained using 
Equation (6.9) and the results of Section 6.3 and are summarized in Table 


6.1. These temperatures together with the electron densities previously 


TABLE 6.1 


AVERAGE ELECTRON TEMPERATURE 


| 
Time (us) Flectron Temperature (eV) | 
' | 
' ! 
| 200 Torr 400 Torr | 600 Torr 
t ay ; i 
| 30 gO OXHe42 | Osdnk 20 Wes Opts tea? | 
| 50 | 0.4 4 1 | O.3itent | OrSetued | 
| 
| 100 | 0.2 + .05 Ovnd++403 oper ae: | 


measured allows one to calculate the parameters Ro a, and y which will 


indicate the conditions under which the probe was operated. This is 


discussed in the following section. 


6.5 Probe Measurements 


6.5.1 Probe Theory 


The high electron density in the discharge implies that the plasma 


sheath around a rotating cylindrical probe is Minin ine this case, cue 


probe current is either due to diffusion or convection of ions into the 


; 2 Dal ; 
sheath. The current is dominated by diffusion if Ro x <1 and by convection 


rR eee In the first case, the probe current per unit probe length 
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pe 1/2. ge LD 


IT=4 
e i . £ 


any a (6.10) 
mere ¢€ is the electron Charge, We the ion mobility, oe the probe radius, 
Ve the probe speed, K is the Boltzmann constant, T, the electron temp- 
erature, and n is the electron density. In the second ease 

poi hesiie 3/4,1/2 


ipa £ UE Rests v8) 


ofoes (6.11) 


where V is the probe voltage and ES is the permittivity of free space. 

The most significant difference between the two cases as seen in Equations 

(6.10) and (6.11) is in the voltage dependence of the probe current. 
Clements, et. isk have measured the probe current in the range where 

a transition from diffusion to convection dominated current is expected. 

In agreement with Equations (6.10) and (6.11), a gradual change from 

2 to v? dependence of the probe current was observed. However, the 


measured currents were up to one order of magnitude smaller than the predicted 


one 
value (Equation (6.10)) when Rio Nose 


Gro rrobe Circuit 
The same probe as described in Section 5.4.2 was used to measure the 


electron density in the discharge. The density could not be measured earlier 


than 50us after firing the discharge because of arcing from the high voltage 


electrode to the probe. The oscilloscope was protected from accidental 


arcs by a metal oxide varistor. Unfortunately, this device had a large 
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capacitance which reduced the bandwidth of the probe circuit. The response 
time was estimated to be 20 us. The probe was biased negatively with a 
battery pack and the voltage was selected with a turn of a potentiometer 
divider. This was bypassed by a 15uF capacitor. The load was varied 

from 102 to 10k depending on the probe current. The probe signal voltage 


was kept much less than the bias voltage; currents were measured as high 


as 90mA. 


6.5.3 Results of Measurement 

A typical oscillogram of the probe current is shown in Figure 6-7. 
As described in Section 5.4.2, the lower trace is a reference signal to 
locate the probe relative to the plasma. Probe current was not very re- 
producible with up to 50% variation from shot to shot. The discharge also 
tended to develop at slightly different locations. Hence, it was necessary 
to take many shots and average the results. Probe current was measured 
from 50us to 1 ms after the discharge was fired. Five shots were taken 
per data point. 

From the known probe speed, the probe current can be recorded as a 


function of position and time. The total current to the probe is: 


ere) I(y) dy 
t fe) 
or 
Fs L(t) rdr 
Peer) = J ————— (6.12) 
ah r -r, 
ah 


where I(r) is the appropriate formula for the probe current per unit 


length. Equation (6.12) is of the same form as Equation (6.5). Hence: 
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which yields n(r). 


BE T, is expressed in degrees Kelvin, and peninglorr;sthensthe 


following relations hold: 


qT. 1/2 
a = 272 (>) 
_ 5.8x10° 
? 
eS 


Peay choad Oxon) 
e 2 
at, 
e 
Sa Toye oe 
Raercex Peal 35x10 ( a ) (6.14) 
Ra 3 p> 
2 Sal, SOLO (eo) 


- cay ? 
The conditions R Novae R, eae: and R, a /y <<l specify the range 
e 


of electron density where Equation (6.11) applies. Equation (6.11) is valid 


when: 
ZY: Ds 
ZOO Torr: 6. 75x10 o<<n<<9. 2x10 /T, 


22, 2 
400 Torr: 8.9x10O<<n<<1.84x10°°/T, 


on 
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Belem 2). (0.0 egktsupl Lepkiqgé (11.8) nokianpa 


4 


600 Torr: 3x10 cen <<2.76x1022/T 2 
e 


From the measured values of n and T.> it is seen that Equation (6.11) 
does not apply up to 100us after the discharge is fired. Furthermore, the 
measured probe currents were about one order of magnitude less than that 
predicted by Equation (6.10). This is in agreement with the results of 
Reference 6. 

For times later than 100us, the electron density and the temperature 
have decayed sufficiently and Equation (6.11) applies. To check that this 
is indeed the case, the probe current was measured at different bias 
voltages. The results are shown in Figures 6-8. The vk 
dependence of the probe current indicates that the current was convection 


dominated. The electron density for t=100us and later are shown in Figure 


6-9 at the three pressures considered. 


6.5.4 Plasma Decay Mechanisms 


The results of the interferometric measurement and the probe measure- 
ment were combined and Figure 6-10 depicts the decay of the electron 
density at r=0. The rate of decay is determined by recombination, diffusion, 


and cooling of the discharge. 


The cooling time due to conduction by neutrals can be expressed using 


vA 
simple kinetic theory as : 


at 


t+ = 10. ECON (6.15) 


where L is the electrode separation, v and A are the mean thermal 
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velocity and mean free path of the neutrals, respectively. tee 300us 
while the experimentally obtained cooling time is approximately 40us (see 
Table 6.1). Therefore, other cooling mechanisms such as line radiation 
operate. Because of the rapid cooling of the plasma, the charge density 
decay is not determined by the cooling rate. 

Referring to Figure 6-10, the plasma decay time is +t = 30us. The 
diffusion time of the ions is approximately equal to is and hence, 
diffusion is ruled out as the dominant loss mechanism. 


. : A : Ae : : ns 
The recombination time is Tp = ewe The recombination coefficient 
R 


varies quite rapidly with density and temperature so it is hard to estimate 


the actual TRe Recombination can account for the observed decay if, 


The above condition yields Op in general agreement with published 


Peles?” Geo seni). Of course, for sufficiently long times, diffusion 


loss will eventually become the dominant loss mechanism. 


6.6 Jamin Interferometry 


A Jamin interferometer was set up to measure the neutral atom 


concentration in the discharge. A He-Cd laser was chosen as the light 


source C0, =441.6nm) . This means that the neutrals as well as the charges 


will contribute to the observed fringe shifts. Since the electron density 


is already known, it is easy to separate the contribution of the electrons 


and ions to the total optical refractivity. The pertinent equations needed 


were already discussed in Section (Ne sip 
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Figure 6-11 depicts the experimental setup while a typical oscillo- 
gram of the observed fringe shifts is shown in Figure) 6-127) Untortunatery. 
it was quite difficult to determine the exact time when the fringes turn 
over (fringe shifts change sign). Therefore only the fringe shifts at r = 0 
were considered reliable because there was no question where the fringes 
should turn over. During the initial part of the current pulse the electron 
density will be increasing while the neutral density decreases. In the 
afterglow, the electron density will be decreasing and the neutral density 
increasing. Both electrons and neutrals will therefore produce fringe 
shifts in the same direction. Furthermore, the contribution of the 
electrons were negligible at the times considered. From the oscillograms, 
the average atom density was computed assuming that the electron, ion and 
neutral gas distributions in the discharge had approximately the same 
boundary. These are shown in Figure 6-13. The uncertainty in the gas 
density is quite large because the total fringe shift for a change in density 
from ambient to vacuum at 400 Torr is only approximately 3.5 fringes (2=lcm). 

There is a factor of two difference between the measured gas density 
and the theoretical equilibrium value at a temperature equal to the measured 
electron temperature. This discrepancy is not surprising considering the 


uncertainties in the neutral density and electron temperature measurements. 
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Figure 6-2 Diagram Used in Abel Transformation 
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Figure 6-3 Mach-Zehnder Interferometer Signal 
The signal from the Ge:Au detector shows typical fringe 
shifts. The sweep speed is 5us/div. The decrease in 


signal is due to plasma absorption. 


97 


7 en 8) oe Le ee ear 
ve 1 | Pr 
: 7! oe ri aa 4 4 ike 
ro 2 oa 
\ i a ws 
re ee y 1 hae 
t 


with smvquonnt aia 
sort a dnatlites weal nosiiast’ ae ‘ie we tomate att 
| 1h Wknadek wit Web a bewegn yore ott sas iia 
| mh) dorsiusi otros ay pat ‘ei purvibe 


rl 
| il ' 
Ie} 
Pl 
hi 
| i} D 
oh i 
We 
| 
i 
r 
| 
\ 


98 


°(9Xeq 90S) eJep ay} OFUT P9IITZ Soet4es Jomod e WOATS ST SUTT PTTOS 24L 


“SoWT], SNoTiea 3e SITUS VBuTAg peAresqg (#) 7-9 2ANSTy 


(wo) JONVLSIG IWIdVs 
Zill_ sO ws 808 WACOM YO. 20 


4401 OOZE=d 


L4lIHS SONIYS 


99 


“souT], snoTze, 3e SIFTYS esuTIg peAresqg (4) 7-9 aan3ty 


(Wo) ZONVLSIG TWwIdvY 
al oF 8'0 90 v0 20 


stool 


440) OOVEd 


L4IHS SONIYS 


100 


“sowT] snoTaeA 3e SIJTUS SBuTA_ peAresqo (9)¥-9 2ANZTY 


(wo) JONVLSIG IIdvY 
Z| ol 80 90 +0 20 


sdool 


4JO0l O09 =d 


LSIHS SONIYS 


gi 

- 
uae 
4 


i? | 
a 
f 
4 
: 

<p 4 


* 


2a 


rae) 


P = 200 Torr 


iS 
© 
x 
ig 

6 

E 

© 
> 
hes 
2 
a KO 
(es 

S 30 
Se) See 
te 
O 
LiJ 
ze) 
lJ 


SSS Eee —EeEeEeEeEEEE——————— EE 


0 O.l Ot)” CCR 04F™ 054 = OGe 207 
RADIAL DISTANCE (cm) 


Figure 6-4(d) Electron Density as a Function of Radius. 
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Figure 6-4(e) Electron Density as a Function of Radius. 
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Figure 6-4(f) Electron Density as a Function of Radius. 
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Figure 6-5 CO 


Absorption by Discharge Plasma. 


2 


ihe upper trace is the discharge current pulse; 
the middle trace is the Ge:Au detector signal showing 
absorption; the lower trace is the detector signal 


when the laser is blocked. Sweep speed is 20 us/div. 
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Figure 6-7 Langmuir Probe Current. The upper trace 
is the probe signal as it passes through the discharge. 
The lower trace locates the probe relative to the plasma. 


Horizontal: 50 ps/div; Vertical: 1 wA/div; p=600 Torr. 
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Figure 6-8 Probe Current Versus Probe Bias Voltage. 


The slope of 1/2 implies that the probe is operated in the 


sheath/convection regime. 
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Figure 6-9(a) Radial Electron Distribution. 
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Figure 6-9(b) Radial Electron Distribution. 
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Figure 6-9(c) Radial Electron Distribution 
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Figure 6-10 Decay of Discharge Flectron Density. 
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Figure 6-12 Jamin Interferometer Fringe Shifts 


d = 0; p = 200 Torr; Sweep speed is 20 us/div. 
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CHAPTER VII 


A CASCADE MODEL OF LASER INDUCED GAS BREAKDOWN 


7.1 Outline 

A cascade model that describes the breakdown of Ar under the in- 
fluence of a cO., laser beam is described. The validity of a classical 
analysis of the laser plasma interaction is discussed. 

The breakdown condition is derived from a square pulse approximation 
taking into consideration the significant loss mechanisms during the plasma 
buildup in the focus. The breakdown condition is assumed to apply to the 
more realistic but complicated case of a time varying pulse. 

An attempt to verify the breakdown condition experimentally is 
described. A charge collector system was used to measure the electron 
density at breakdown. The results were inconclusive, however. 


Numerical solutions of the cascade process are performed. The agree- 


ment with experimental results is good. 


#22 Introduction 
The problem of cascade breakdown under the influence of light was 


first theoretically investigated by Zel'dovich and Raizer . In their 


fundamental paper, the quantum kinetic equation for the electron distribution 


. s 2 Ww . 
finction was solved in the "classical limit". The energy absorption 


mechanism was inverse bremsstrahlung absorption by electrons during collisions 


with neutrals. In the classical limit the rate of energy absorption via 


inverse bremsstrahlung is the same as that used in microwave theory, namely: 
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The validity of the classical analysis was investigated by Friedland” 


He showed through particle simulation that for CO, laser breakdown of Ar, 
classical theory is applicable even though the conditions set forth by 
Zel'dovich and Raizer on the constancy of the collision frequency is not 
met. Consequently, the breakdown of Ar by an intense co, 
by the classical Boltzmann equation with the appropriate collision terms 
which will be discussed later in this chapter. 

Aside from Reference 1, a good presentation of laser cascade break- 
down including detailed derivation of absorption coefficients can be found 
in Reference 3. 

Assuming there are a few initial electrons in the focal volume, the 


build up of electrons may be described by a differential equation for the 


electron density, n: 


C701) 


where v is the electron gain rate due to energy absorption by electrons 


via inverse bremsstrahlung and Vo is the loss’ rate due to electron 


diffusion out of the focal volume and energy losses during collisions. 
The solution of Equation (7.1) gives the electron density as a function of 


time n=n(t). If breakdown occurs at the end of the pulse, t= with 


n=n then ny =n(t) The breakdown threshold, therefore, is the laser power 
D- : 


i ive at the end of the pulse. 
Pax reduired to gi ny 
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In order to solve Equation (7.1) several things must be taken into 
consideration. First, the density at breakdown ny is unknown. Secondly, 
the loss and gain rates are not constant. Hence, as the electron temperature 
varies with the instantaneous laser power, the rates will be changing in a 
complicated manner. The rates are also dependent on the electron density. 
Thirdly, the plasma might not be in equilibrium, i.e., the distribution 


function is non-Maxwellian. This will affect the gain and loss rates. These 


are considered in the following sections. 


7.3 Breakdown Condition 

Although the density at breakdown, n, is unknown it is still possible 
to obtain a condition for ny, at breakdown. There are two possible 
energy gain mechanisms: electron-ion inverse bremsstrahlung and electron- 


neutral inverse bremsstrahlung. During the early part of the charge build 


up, the electron density is low. For this reason the electron-ion inverse 
bremsstrahlung term was neglected in previous work . In References 4-8, 


breakdown was defined arbitrarily as full ionization or in some cases, 


a =—5 4.5 ; 
n,=10'°cm™ with an initial density, a =10"cn . Stiil others have 


defined breakdown as 43 electron generations in the cascade It will be 


shown that the electron-ion term is important in that it sets a condition 


for an effective value of n° It will also be shown that although the 


electron-ion term is important, it can be disregarded when solving Equation 
(7.1) for the charge build up thereby simplifying the differential equation 


as long as the breakdown condition to be derived is used. 


Define v as the gain rate due to electron-neutral inverse 
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Here, US is the ionization potential of the gas, E is the rms electric 


field strength of the laser, e and m _ the electron charge and mass, 
respectively, ee is the collision frequency for momentum transfer and w 


is the laser angular frequency. The absorption coefficient corresponding 


to this gain rate is: 


Here, F is the light flux. The electron+ion inverse brem- 


ee a 
sstrahlung absorption coefficient is : 
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Equation (7.6) reduces to the case of no electron-ion inverse bremsstrahlung 
absorption when the density is low so that the second term on the right 


hand side of the equation becomes much smaller than the first. 


Assuming further that E, ye and Vv, are constants, the solution of 


Equation (7.6) is: 


n_exp(v tay, )t 


n(t) =— 7 (7) 
=n) afexp(v, —v ,)t-1} 
where 
petox00 oe 
a = 
1 
oF “Vv, 
The density at the end of the pulse is: 
n_exp(v res Ms 
n(t) = MOMMA A ee Cis) 
1 { (v — \t-1} 
-nalexp yeny 
Taking the logarithm of Equation C758): 
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Figure 7-1 is a plot of (varv9)t versus n(t) with n=10°cm> and a=0, 
Wa7x10"’. and Teoexi0" For a given laser power and gas pressure, the 
density at the end of the pulse, n(t), will lie along the curves shown in 
Paguce /-1, e.¢g., point A on line AB. If the laser power is increased, 
n(t) may correspond to point B. If electron-ion inverse bremsstrahlung 
were disregarded, the solution will be along the straight line AD. Note 
that point C corresponds to the same (v =v9)t as point B and in fact 
approximately corresponds to any point ners n(t)>l/a. Therefore, the 
electron-ion term can be disregarded as in previous work if the calculat- 
ion is terminated at n(t) = l/a = n- Returning to Equation (7.1), its 
solution assuming constant Ev and v, is: 


£ 


Ai) = ngexpt (v,"~¥9)T) (7.10) 


Setting nan (t) = 1/a in Equation (7.10) is seen to be approximately 
equivalent to setting the denominator of Equation (7.8) to zero: 


Lie vk Be at) 
1 - an texp(v, V,)t Ly =20 


This is a sufficient condition for breakdown. 


The above discussion leads to the following conclusion. When 
modelling the buildup of electrons in the focus during breakdown, the 
simplified problem of no elect 
to the more realistic case which takes it into account if the density at 
This is true only if the actual density 


breakdown n is taken to be 1/a. 


at breakdown is larger than 1/a. 


ron-ion inverse bremsstrahlung is equivalent 
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i 
vo =v 
The factor 2+ =~ 1. Hence, 
v 
g 
fe 3 43x10 (7.12) 
b ; Cc ; 


For argon with se 2-3eV. ee seals where p is the gas pressure in 


Bocr.” Therefore: 


mn, = 9.4x10"> p (7.13) 


dite oe : 
It can be seen that for moderate pressures, ni < 10 “cm ~ which was used 


by others in previous calculations. The breakdown condition is pressure 


dependent and is equivalent to prescribing the degree of ionization of the 


gas: 


ar 


12 


0.277 


= 


gas 


7.4 Charge Collection Results 
In an attempt to experimentally determine ny» a charge collector 


consisting of two flat parallel plates 5cm by 5cm was constructed. The 


laser was focused midway between the two plates. Figure 7-2 shows the 


equivalent input circuit seen by the collector. It can be easily shown 


following the analysis of Reecher ss: that the output voltage signal is: 
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where 
R_R 
iva 
R_ = a> R 
2) RR, 


I(t") = collector current 


C = collector capacitance 
Cc, = input capacitance of electronics 
C+C = 100 pF 


If the collector field is high, Townsend ionization will take place 
as the electrons drift towards the anode. Let a be the first Townsend 
coefficient, No the number of electron-ion pairs created by the laser, 

d the distance between collector plates, v_ and Vy be the drift velocities 
of the electrons and ions, respectively and T_ = SUAS? oe TL = d/v,. Since 
the charges are created midway between the plates the currents due to 
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5 = exp[-52/(E/p)] (7.18) 
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It can be seen that the signal should increase rapidly with E/p due to 
the increase in a/p. 

At very low E/p, the collection time is so long that losses can 
take place. The electron diffusion time is approximately: 


t 2A? /D_ 


d 
where A is the diffusion length and D is the electron diffusion coeffticienc, 
Taking the ratio of i) and £ 


ada 


7 q2 ak 

=o a C25) 
eat 

where T, is the electron temperature in eV and V is the voltage between 


the collector plates. Taking d = 1.5cm A = L/w = 1.6 cm, then: 


dk a3 
| See 


Hence, except for very low V, diffusion loss can be neglected. 

Figure 7-3 (a) and (b) depict the collector characteristics for 
p = 30 Torr and p = 50 Torr and laser energies E = 2.5J and 1./5J, resp- 
ectively. For E/p less than 4 V/cm-Torr, the fieldvcannot collectsals 
charges while for E/p greater than 15 V/cm-Torr, Townsend ionization is 
dominant. A collector field of 10 V/cm Torr was selected as the operating 


point. 


The number of electron-ion pairs created by the laser at threshold 


were measured for both cases of when a spark was produced and when a spark 
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was not produced. When a spark was produced, the number of collected ions 
was consistently 10°? regardless of the pressure, as shown in Figure 7-4 (a). 
When a spark was not produced, the number collected varied considerably 

from shot to shot. This was also observed by Tulip and Sept inn at The 
smallest number of electron-ion pairs collected was tee as shown in 

Figure 7-4 (b). This is an order of magnitude less than that obtained by 


Tulip and Seguin in air. Obviously, the true density at breakdown should 


be between 10° and 102 With a focal volume of So. cme these limits 


correspond to ates and oa ene: 


The breakdown condition discussed in Section 7.3 is valid only when 


3 


1 - Uo ge fea 
the actual density at breakdown (between 10 Hes and 10 cm §) is greater 


than 1/a. To check whether this is satisfied, assume n= 


ne | i =9 ; 
oe om e The build up from Peele TO? mee ak6) 3 cm needs approximately 20 


pe eel Dears 
cascade generations. If the gas is preionized such that ny =10° “cm , only 
4 cascade generations will be required to break down the gas. Since the 


required laser power to break down the gas is approximately proportional 


to the number of cascade generations required, a threshold lowering of 


4 ? ; 13 -3 
= 20 7 0.2 should be observed if the density at breakdown is 10 cm . 
18 


On the other hand, if Dee 10 cee 30 cascade generations are 


5 -3 i f = LOWS 
needed from Ns 10°cm ~ while it takes 10 cascade generations trom fei 


= 10 ; 
cm 4 Hence, a threshold lowering of =~ 40.0 0.33 should be expected in 


this case. 


Measurements of the breakdown thresholds reported in the next chapter 


12. -3 ? ; 
indicate that with a preionization of = 10°“cm ~, there is a lowering of 


the threshold to only 0.4 its value without preionization. This indicates 


Loa h 10l3em > 
that the actual density at breakdown is closer to 10 cm than . 
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Hence, the condition discussed in Section 7.3 is satisfied. The 
breakdown condition obtained in Section 7.3 will be used in the following 


sections. 


7.5 Boltzmann Equation 


As previously discussed, the pulse shape, cross sections and varying 
temperature make ae and DF, variable in time. Expressions for ve and 
Vv, as functions of time (laser power) and electron temperature are needed 
before Equation (7.1) can be solved. This means that a rate equation for 
T, is also needed. These are obtained from the Boltzmann Kinetic Equation. 


The Boltzmann equation that describes the time evolution of the elect- 


i 
ron distribution function has been derived by Allis : 
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engor e  ge ge | Opie 1/2 ) 
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: 7.26 
-) Ay (e)f(e) - A, (e)f Ce) (7.26) 
alt 


where 


£ = distribution function 


m = electron mass 

M = ion mass or atom mass 

Vv = collision frequency for momentum transfer 
c 


A = excitation rate of level i 
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A= 1onization rate 


Baten hig ee sntaual! iy nt) » ih i tone 


& 
*.. 


EL ‘wisn : 


iat ntiawohe whe kee sc whit 
4 wat Geto) Ba nga ai ih a a hs Pou nt saith te 


Mae i ee ey pat be | i 


129 


Q 


e 
= heating rate = — 
Jj & m 


The first term gives the changes in f due to diffusion. It is the only 


space dependent term and so the following approximation was taken: 
2 
Vf =- £/A (ieZ7) 


Also, define: 


2€ 

2 

3mv (e)A 
c 


Ay Ce) = 


The second term is the heating rate and the third term is the elastic 
collision term. The last two terms give the losses due to excitation and 
ionization. The term under the integral sign is due to the appearance 
of electrons of energy e¢ after exciting or ionizing collisions of electrons 
of energy €'>c. Recombination losses are neglected as well as hydrodynamic 
effects, photoionization, radiation loss, collisions of the second kind, 


ionization of excited atoms by electron collisions, and heat conduction. 


ioe Integral Term 


7 , 
Following the analysis of Chan and Moody , consider first the rate 


! 
of appearance of electrons of energy € due to electrons of energy € 


te ‘ ; € 2 BYyrthe 
exciting an atom during a collision, (excitation energy x,? y 


conservation of energy: 
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Hence, 


Ge ve!) fe") | = y Be ace JA, Ce" ie") C729) 
iit i 
ex 


Similarly, electrons of energy e«' can ionize an atom creating an electron 


of energy ey and another of energy Ey where: 


= q pe = e ° ° ° 

ey + Ey iS er (e, ionization potential) C/230) 
Bet P(e, e'-€,) be the probability that es- € and E, = e' - €, -€ and 
P(e' - er 7 e,e' - Ey) be the probability that E,= € and Ee Pe e' — Ee, — &- 
me, c' = Ey) satisfies the normalization condition 

e'-€, 
f Rc.ee = €,)de =] O73 1) 
oO 


and 


P(€,e'-€,)=0 if e<O and e>e'- Ey 


The lower limit in Equation (7.31) corresponds to the case when e.= 0 


and the upper limit to the case when E5= Ot, Erhen: 


Ge...) f Ge") = A, (e")£(e') [P(e,€'-€7) Be P(e'-e,-€,€'-€,)] 
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Adding the excitation and ionization rates, the total rate is: 


eee ite) = £(e')[2A,(e")D(e,e'-e,)+ ) A ke eet ae) 
ab au 


Ch533) 


lL 
where D(e,€'-€,) = alP(e,e! - Ey) +P Ce = eu Clee ee) 


it 2 


; . un 
which is symmetic about ¢€ = ae! - Ey) because each electron of energy 
€ produced during an ionizing collision must be accompanied by another 


of energy e' - e;- D satisfies : 


a Dte,e' = ede =] (7.34) 


7.5.2 Conservation Laws 
It can be shown that the excitation and ionization terms discussed 
in the preceding section satisfy the conservation laws. Integrating 


meess)f(e') over all € and. c', 
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Changing dummy variables from c' to e, 


f Gle,e")£(e")de'de = S£(e) [2A (e) + ) A (e) }de (7.36) 
a at 


This is the law of conservation of mass for exciting and ionizing collisions 


involving electrons of energy e'>c. Similarly: 


ete t(¢ ')detde = f(e-e )E(e)A (e)de + ) fA. (e)(e-e, )E(e)de 
af i i 


Kot) 
which is the corresponding law of conservation of energy. 
7.5.3 Density and Temperature Equations 
The relations: 
n= / f(e)de C7238) 


and 


eer = tee (ede 
2 e 


are used to obtain the rate equations for the electron density, n and 


electron temperature, T.. Integration of Equation (7.26) yields: 
e 
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at Apfedt(e)de + sf A, (e)£(e)de = (A, - A))n (7.39) 
p) 5 
pets uh) = - f eA (e)t(e\de = 7 2 23/2 ce esl) 
2s Bee ede ee Os ean gee 
€ 
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- ) e, SA, (e)£ (ede 
i a 
Or 
oT 
pee ey hk A(T) 
dt 3 D J el teed 2 e 
= 3 as 
Re a4 140 
ey i 2 ae Ay] 
where ns ADs and so forth have been defined such that for example: 


nes 
A, os oF A, (e)f (e)de 


— 


A, is the average ionization rate, while Ay is the average diffusion 


rate, Ry is the joule heating rate, aS the elastic collision loss rate, 


— 


Ry the average energy loss due to diffusion, ae the energy loss due to 


ionization and c¢ A. is the energy loss due to excitation. Equations 
a gle 


(7.39) and (7.40) are the desired rate equations. Its solution will give 
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n=n(t) and the threshold will just be the smallest laser power that will 


give n=n, . 


7.5.4 Maxwellian Distributions 
Because of the deep Ramsauer minimum in the collision cross section 

of argon, it is expected that the distribution function is not Maxwellian 

but one which gives greater weight on the low energy side. A numerical 

solution of the Boltzmann ecnet ion and particle imitans one show that 

this is indeed the case. However, the distribution function was not 

very different from a Mexweiiien. Therefore, the terms on the right hand 

side of Equations (7.39) and (7.4) derived from a Maxwellian distribution 

are not expected to be very different from their true values. Furthermore, 

the greatest contribution to these terms come from electrons in the tail 

of the distribution function. The distribution functions obtained in 

Reference 3 are Maxwellian in the tail. 
Assuming a Maxwellian distribution: 

1/2 -e/T, 


f(e) = 2n( — >) e 
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(7.41) 


the collision rates become (energy in eV): 
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where P is the laser power per unit area in the focus. 


The collision frequency is defined as: 


where o is the cross section and Nes is the electron velocity. 
m 


Equivalently, On is: 


v = 1.92x10°" ell? o (e)p C7243) 
c m 
Where c¢ is in eV, oF in aie and p in Torr. Similarly, 
2 7.44 
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Here, oF and me are the electron impact ionization cross-section 


and the excitation cross-section, respectively. The appropriate 
expressions and the approximations used in the calculation including 
reference sources are summarized in Table 7.1. The rates are plotted in 
Figure 7-5 (a) - (c). A:polynomial up to an 8th degree was fitted to 

the rates and used in the solution of the density and temperature equations 
(Equations (7.39) and (7.40)). 

As the electron density builds up, a point is reached when the 
diffusion loss will be greatly reduced because of ambipolar effects. To 
take this into consideration, the diffusion loss terms Ay and m were 
multiplied by a function F which takes on values from 1 when n is 
low to D,/D. when n is large. The functional form of F was similar 


to that prescribed by Allis and Rosena ae 


for a plasma between two 
Z ; : : 
parallel plates. F plotted as a function of nA iT. is shown in Figure 


H-6. 


7.5.5 Numerical Solution 


A finite difference method was used to solve Equations (7.39) and 


(7.40). The following substitution was used: 


Jtlen J 
Bie 1.46 4 (7.46) 
ot AC 
on diy attend (7.47) 
at at 


where pitt is the present time. The terms on the right hand side of 


j “ 
Equation (7.40) were first evaluated at the time t-. Equation (7.40) 
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was then used to evaluate and AL and Equation (7.39) yields nitt, 


: j 
+1 n 
n? = e+ (7.48) 
| oa At 
where et = A, = Mais The procedure was then repeated using the 


new values of T, and n. Two or three iterations were found to be 
sufficient for convergence. Accuracy was up to three significant figures. 
The time step size selected was At = 5100s: second. Other schemes 
such as the 4th order Runge-Kutta and the Adams predictor-corrector 
method were tried. No significant advantage of these methods over the 
one used was found. A program listing is given in the Appendix. 
The laser pulse shape was approximated by that shown in Figure (7-7). 
This should be compared with the actual pulse shape given in Figure 4-2(b). 
The procedure described in Chapter III insured that priming elect- 
rons are present in the focal volume just before the laser pulse. Since 


3 


= ree 
the focal volume is ~10 Dear an initial value of n=10° cm ~™ was chosen, 


i.e. there is one electron in the focal volume to start the avalanche. 

The laser peak energy was adjusted and the difference equation 
solved. The value of the diffusion length is still unknown and how this 
is set is discussed in the next paragraph. If for a given laser peak 
power the gas breaks down, the power is reduced by 10% and the differential 
equation is solved again until there is no breakdown. The threshold is 
taken as the average of the smallest power with breakdown and that without 
breakdown. 

To determine the value of the diffusion length, A, the breakdown 


threshold was measured using the procedure described in’ Chapter TIT, 


From 50 Torr to 10,000 Torr. 
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assuming Gaussian beam optics was used to calculate the experimental 

peak power per unit area at threshold. A was then adjusted until 

there was reasonable agreement between the experimentally measured 
threshold without preionization at p=100 Torr and the results of the 
numerical calculation. The best fit was obtained with A=67um. This is 
larger than the theoretical value of 100um/2.405=41.6um and this implies 
that the plasma is not confined to the focal volume. The experimental 
thresholds without preionization and the results of the numerical modell- 
ing are shown in Figure 7-8. 

The time breakdown occurs depends on the pressure. At low pressures, 
breakdown occurred very near the peak of the laser pulse. At very high 
pressures, on the other hand, it occurred very near the end of the pulse. 
This was due to the very high particle loss rates when the pressure was 
low. The time of breakdown plotted against the pressure at threshold is 
shown in Figure (7-9). 

The peak electron temperature varied from 2 to 4 eV. Higher electron 
temperatures were obtained at low pressures, again due to the nature of 
the losses. This is also shown in Figure (7-9). Figure (7-10) depicts 
the laser pulse, electron temperature, and electron density as a function 
of time when breakdown occurs and when it does not,at p=100 Torr. The 
electrons are heated rapidly after passing through the Ramsauer minimum. 


Thereafter, the electron temperature follows the laser pulse. Note also 


that a small change in the laser power leads to a large change in the 


peak density generated in the focal volume. This is an indication of the 


sharpness of the threshold. For comparison, the electron density 7burcup 


at 6000 Torr is also shown in Figure (7-10). 
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Since the electron and neutral densities under experimental 
conditions are known, the breakdown thresholds can be measured and compared 
with numerical calculations. If reasonable agreement is obtained, not 
only has the cascade nature of breakdown been confirmed, but more importantly, 
the possibility of using breakdown as a diagnostic, especially for neutral 


atom density measurements , appears very promising. 
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Figure 7-3(a) Charge Collector Characteristic 
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Figure 7-3(b) Charge Collector Characteristic 
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Figure 7-4(a) Number of Electron-Ion Pairs 
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Figure 7-4(b) Number of Electron-Ion Pairs 


Collected when a Spark is not Formed 
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Figure 7-5(a) Laser Heating Rate as a Function of Electron 
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Figure 7-7 Laser Pulse Shape Used in Numerical Calculations 
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CHAPTER VIII 
BREAKDOWN THRESHOLD MEASUREMENTS: COMPARISON WITH 


CASCADE THEORY PREDICTIONS 


8.1 Outline 

The co, laser induced breakdown threshold of Ar in the presence 
of an electric discharge was measured with the time delay between the 
initiation of the discharge preionizer and the laser pulse being varied. 
The co. laser was focused inside and outside the discharge. Outside the 
discharge the threshold was lowered due to preionization in the focal 
volume by uv photons and compression of the gas during the passage. of a 
weak shock wave. The measured thresholds were related to the previously 
determined electron density outside the discharge. It is shown that 
cascade theory adequately accounts for the experimental thresholds. 

The lowering of the breakdown threshold during the passage of a 
compression wave through the focal volume is observed. Neutral atom 
densities corresponding to the amount of lowering are compared with pressure 
probe measurements of the shock strength. 

When the laser was focused into the discharge region the changing 


electron and neutral gas densities combine to alter the breakdown threshold. 
The possibility of using breakdown to predict either the electron or the 


neutral density is explored. 


oa2 Experimental Setup 


The experimental setup has already been described in Chapter IV. 


Figure 4-1 is reproduced as Figure 8-1. A 5cm NaC& lens was used to focus 


e 
the laser in the breakdown chamber. Since the probe measurements were don 


rt, a 
Mithla metal plate shielding the discharge except for a small uv po 
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similar plate was placed in the breakdown chamber when measuring the 
thresholds outside the discharge. The laser was focused into the center 

of the approximately cylindrical plasma. The laser was pulsed at a rate 

of one pulse per second as in the breakdown measurements without deliberate 
preionization. This was done more for convenience than necessity since the 
procedure described in Chapter III is not necessary when there is 


appreciable preionization. 


8.3 Breakdown Measurements Outside the Discharge 


The breakdown thresholds in Ar were measured at the same pressures 
and distances from the discharge as the probe measurements were made. The 
results are presented in Figures 8-2 (a) - (d). The thresholds have been 
normalized by their values without preionization. The initial time, t 4-0 
was arbitrarily set to be 4us after the initiation of the discharge pulse. 

The lowering in the thresholds was due to the preionization of the 
gas by uv photons emitted by the discharge. As the photoplasma decays a 
gradual increase in the threshold was observed. Aside from the lowering 
due to the preionization, the passage of a shock wave in the focal volume 
was detected and manifested as a further lowering of the breakdown threshold. 
These can be seen in Figures 8-2 (a) - (d), and will be discussed later 
in the chapter. 


Excluding the areas in the measured thresholds where shock effects 


were present, the experimental values of the thresholds were plotted against 


the initial electron-ion density in the focus. This is shown in Figure 


8-3. At the same time the numerical model discussed in Chapter VII was 
id 
used to generate theoretical threshold values and are plotted as soli 


h eS 
lines in Figure 8-3. The experimental points do not all ‘fall on the theore 
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ical curve. The discrepancy between the two may be due to a higher 
initial density in the focus in the case of no preionization than pre- 
viously a osuned.. Uncertainties in the electron density, spot size, 
experimental thresholds, and other factors should also be taken into 
consideration. 

The bend in the curve especially at low pressures is due to a 
transition from free diffusion loss to ambipolar diffusion loss of electrons 
from the focal volume. This transition was first observed by Brown and 
erith: . 

The weak dependence of the breakdown threshold on the initial 
electron density as indicated in Figure 8-3 and Figure 8-5 makes it very 
difficult to obtain accurate measurements of the electron density from 
threshold measurements. A small error in the measured threshold will 
correspond to a large error in the density. On the other hand, the thres- 
hold varies quite rapidly with the neutral atom density (see Figure 7-8) 
except at very high pressures (10,000 Torr). Therefore, if the electron 


density can be measured independently, neutral atom concentrations can be 


obtained with reasonable accuracy. Errors in the electron density will not 


be very critical because of the weak dependence of the threshold on it. 
With the plate shielding the discharge removed, the threshold was re- 
measured at the time of arrival of the expanding shock wave at the focus. 


The plate was removed in order to obtain measurements very close to the 


discharge. The results are shown in Figures 8-4 (a) - (£). The thresholds 


e 2 ik 
without the plate are in general lower than those obtained with the meta 


shield. Furthermore, it takes a longer time for the plasma to decay. 


These indicate a higher degree of preionization without the shield. The 
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higher preionization may be due to more efficient use of the photons 
emitted by the discharge. The effective photons may not be emitted 
uniformly from all regions in the discharge but mainly come from particular 
areas such as near the electrodes. Photons from these areas may have been 
partially blocked by the plate. There are no measurements of the electron 
density due to uv radiation for this case because the probe cannot be 
exposed to uv. Therefore, it was assumed that the electron density could 
be obtained with sufficient accuracy from the measured thresholds. From 
the known gas pressure and breakdown threshold just before the shock, the 
electron density can be read from Figure 8-5 which gives the theoretically 
calculated thresholds as a function of the initial electron density and 
pressure. The electron density during the passage of the shock is then 
assumed to be the same as its value just before the shock. With this initial 
electron density and the new threshold in the shock, the gas pressure and 
hence, the neutral density can be obtained from Figure 8-5. Interpolation 
between two curves may be necessary. 

The pressure jump in the shock has been measured previously (Chapter 
V). Assuming the equation of state of the gas is known (perfect gas) and 


assuming further that y = 5/3, the gas density jump in the shock can be 
computed using shock theory. The results are the solid lines in Figures 
8-6 (a) - (c). The data points are from the breakdown measurements. It can 


be seen that for this case, breakdown threshold is a direct measure of the 


neutral atom concentration of the gas. 


8.4 Breakdown Measurements Inside the Discharge 


The breakdown threshold when the laser was focused in the center of 
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the discharge was measured as a function of time delay between the dis- 
charge and laser pulses and the as pressure. The threshold ratio can 
exceed Breas This occurs when there is a reduction in the neutral density 
which in this case is due to heating of the gas. As the gas cools down, 
the threshold slowly returns to its initial value. This can be seen in 
Figure 8-7 which depicts the measured thresholds at the center of the 
discharge. This can also be seen in Figure 8-4 (a) at approximately 200us 
after the discharge was fired. At this time the discharge has expanded 
to the region of the focus. Although there is also an increase in the elect- 
ron density as the discharge expands into the focal volume, the effect of 
the reduced neutral density predominates because of the weak dependence of 
the threshold on the initial electron density. 

From the measurements and the known electron density (Chapter IV) the 
neutral atom concentration in the discharge can be computed and compared 
with the results of the interferometric measurements. This is shown in 
Figures 8-8 (a) - (c). Also shown in the figure are the expected neutral 
densities if equilibrium is assumed with a temperature equal to Ty obtained 
from the absorption experiment. 

It can be seen from Figure 8-8 that the computed neutral densities 
from the breakdown measurements are consistently smaller than those obtained 
from the interferometric or absorption measurements. This may be due to 
the difficulty in visually observing the laser spark superimposed on the 


bright arc discharge. In fact, only the tail of the spark which comes out 


of the discharge is observable. It is then very probable that the thresholds 


obtained were larger than the actual thresholds; more energy is required 


to produce a spark which can be observed amidst the bright background. This 


sets-a limitation to the technique which could possibly be overcome if the 
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threshold were defined in a different manner. The most likely candidate 


for this would be a measurement of the transmitted laser pulse. 


8.5 Discussion 


The preceding sections indicate that cascade breakdown can be used 
to measure electron and neutral atoms densities in rapidly changing 
plasmas. It is more effective in cold, low density plasmas (e.g. photo- 
plasmas) where radiation does not interfere with the observation of the 
spark. 

If breakdown is used to measure electron densities (provided the 
neutral density is known), the uncertainty in the results will be large 
because of the very weak dependence of the threshold on the initial density. 
Therefore, breakdown possibly gives only order of magnitude estimates. 

If the electron density is known, rapid changes in the neutral gas 
density can be monitored by observing changes in the gas breakdown threshold. 
Since $v =o a the threshold is more sensitive to changes in n than 
n. The results will be in close agreement with measurements using other 
well known techniques as long as there is no difficulty in observing the 


laser spark. The uncertainty in n, obtained will approximately be of 


ig 
the order of the uncertainty in ¢ since ¢ V7. 
n 


Aside from the limitation of the technique to non-luminous plasmas, 
the magnitudes of the electron and neutral density will impose limitations 


on the practicability of using breakdown as a diagnostic. For very high 


au ld : 
initial electron densities (= soem ), it would be very difficult to 


control the incident laser energy to obtain probabilities of breakdown between 


0% and 100%. An absolute limit to n is, of course, n,- Furthermore, 


ities 
since avalanche breakdown is assumed to take place, only neutral densit 


ay i if 


eee hk 


104 CTS OVER oles sion iia 


‘~~ 
1 


r : ; ; F My . 
PRAY ene Wee Me saat aM as di be ’ 
Ae i} A a Ree iat ' y) i ae i 


Sa be bo es wa Hi ia Ah wot es He ; 


ana m ty) yy Lay Bae ve agi! ein + ie sto M REN tt . om a 


edety Leet La eat oe eee ee nb 


pad Veh 8 Dale ae ae oes 
Des dawn tht oh’ Pe a ae iia ae ! oh stl : 
hes eK, yaaa me PPR Ahi &, a Ait 5 sl tye i Bi) “a 


anehy ot pubes eget ay a Lie. ‘sie i ‘i 


a i fd fi 
&! 4 ‘ed Af) toa & i f : » a ; ; h a! Rr: a4 
ian Dato fs ve ‘ \ Mig ‘ q Baal! 
r 4 a re De eae a it gia we ape 
hy aie ; 


» be “et adiintienrtdere Gidbw Getic thas ae aes laa one: 
; 1% ‘ f ; ive on yy 4 i . i 


} - 


abate ee cam a 194 Hera 


C4 
= 


’ 
; - 
« Be ich me fi Hey i’ 4 t f  d 
7) 
f me 
i 
tal bp} ; oRirag Pil a) aRAT “aad te 
4 i | 
7 va am ’ 
« 7 or ie ie 
N eli ALANS 


161 


£ 17) =3 
greater than ~ 3x10 cm (p = 10 Torr) can be measured. At very low 


pressures avalanche breakdown does not occur. 

As an example of the use of breakdown to determine electron and 
neutral densities, the breakdown threshold of Ar at a pressure of 600 Torr 
was measured 8 mm from the center of the discharge. This is shown in 
Figure 8-9. Since the measurement was performed at a distance much closer 
than those reported in Figure 8-4, effects of the discharge expanding into 
the focal volume should be observed. Measurements of the threshold at other 
pressures and at slightly different distances give similar results. The 
events that take place in the focus are as follows: For very early times 
the discharge emission photoionizes the gas which will be detected as a 
lowering in the threshold. The discharge creates an expanding shock which 
passes through the focal volume approximately two microseconds after the 
end of the current pulse. This is observed as a further lowering in the 
threshold. Meanwhile, the discharge itself is slowly expanding. At approx- 


imately t, = 150us the hot discharge has expanded up to the focus. The 


d 


higher temperature means a slightly lower gas density compared to the back- 
ground density and a slight increase in the threshold is observed. The 


electron density will also increase at this time but the effect of rare- 


faction is more pronounced. As the gas cools, the electron density decreases 


and the gas density increases. The threshold gradually approaches its value 


without the discharge. 


The electron and gas densities determined from the thresholds shown 


in Figure 8-9 are depicted in Figures 8-10 (a) - (b). From t, = 150 to 


400us the electron densities were obtained from the probe measurements. The 


initial decay of the plasma is recombination-controlled with a recombination 


coefficient of ReNOe comes The electron and neutral atom densities 
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follow the qualitative picture discussed above. 

Finally, the photoplasma density at t 4-0 obtained from the breakdown 
measurements was plotted versus distance from the discharge as shown in 
Figure 8-11. In this case, the neutral density will be equal to the 
ambient gas density. As in the previous probe measurements the variation 
of the initial density with distance is not incompatible with an ne 


dependence as shown by the solid lines in the figure. 
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Figure 8-5 Computed Laser Peak Power at Threshold as a 


Function of the Initial Electron Density 
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Neutral Density in the Discharge Versus Time 
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Figure 8-8(b) Neutral Density in the Discharge Versus Time 
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Figure 8-8(c) Neutral Density in the Discharge Versus Time 
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Figure 8-10(b) Variation of Neutral Density with Time 8 mm 


from the Center of the Discharge 
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CHAPTER IX 


SUMMARY AND CONCLUSIONS 


The effects of changing initial values of the electron and neutral 
densities on the co, laser induced breakdown of argon has been studied 
experimentally and theoretically. In particular, the thresholds were 
measured when the laser was fired in conjunction with an electric discharge 
Similar to those used in uv preionized lasers. The experimental setup 
consisted of a co, laser-amplifier system, a breakdown chamber, and the 
electric discharge preionizer. 

The characteristics of the laser system were presented. As well, 
the electrical characteristics of the discharge were studied and the energy 
in the discharge converted to mechanical energy was measured. 

The level of preionization due to uv photons from the discharge 
was measured with a rotating Langmuir probe. To interpret the experimental 
data, present theory on convection currents to probes as published in the 
literature was extended. The radial dependence of the measured photoplasma 
density was not incompatible with an ga dependence which seems to indi- 
cate a two-photon ionization mechanism. This problem is beyond the scope 
of the present work and points to the need to study the photoionization 
process taking place in greater detail. 


The electron and neutral densities in the discharge were measured 


i laser and a He-Cd 
using interferometric and probe techniques. A cO., a 


i dons... fi 
laser were used as light sources in the interferometric investigat y 


output, the fringe shifts 
virtue of the long wavelength of the co, laser put, 


LELeSs. 
observed in a Mach-Zehnder interferometer directly yielded electron densi 
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The probe was used to determine the electron density later in time. The 


decay of the plasma was recombination controlled. 


An absorption experiment was set up to determine the electron temp- 
erature. The plasma was found to codl much more rapidly than would be 
expected from heat conduction losses. This indicated that other heat loss 
processes were occurring, most probably line radiation. 

A cascade model of the breakdown process was developed taking into 
account the laser pulse shape, diffusion losses, elastic and inelastic 
collision losses, and variable collision frequency. Numerical calculations 
were performed based on this model and the breakdown thresholds were 
calculated as a function of the initial electron density and the gas 
pressure. The large lowering in the thresholds in the presence of a dis- 
charge was seen to be due to sufficiently high uv preionization such that 
the main loss mechanism in the breakdown process was ambipolar diffusion. 


The transition from free to ambipolar diffusion loss typically occurred 


around n = 10' om >. 


The breakdown thresholds were determined both inside and outside the 
discharge. Thresholds were also measured without preionization. In fact 


the results of the measurements without preionization were used to check 


the validity of the cascade model discussed above. The gas was first con- 


ditioned to ensure that the measured thresholds were true avalanche thresholds. 


A comparison of the measured thresholds with preionization and shock comp- 


i nerall 
ression with theoretical predictions was made and there was ge y 


laser 
reasonable agreement. This further supports the cascade theory of lase 
induced gas breakdown. 


When the laser was focused in the discharge, a discrepancy between 


i attributed 
theoretical and experimental thresholds was observed. This was 
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to difficulty in experimentally determining the occurrence of breakdown 
when there is strong light emission from the plasma being studied. If 

other ways of determining breakdown can be used this difficulty may be 

overcome. 

The possibility of using breakdown to measure the electron or neutral 
density of a plasma was explored. The threshold was not very sensitive to 
changes in the initial electron density which implies that the uncertainty 
in the densities derived from the threshold is large. It appears to be 
applicable only for order of magnitude estimates. 

On the other hand, as long as n is known and there is no difficulty 
in observing the spark, the breakdown threshold yields neutral densities 
in good agreement with well established techniques. There seems to be 
enough reason for further study of the breakdown process. As a probe it 
will offer a high degree of space and time resolution and the convenience 
of minimum optics and very little or no detecting instruments. The dis- 
advantage seems to be the need for many "shots" before the threshold is 
determined, the dependence of the threshold on numerous parameters such as 


lens focal length, laser divergence and most of all the dependence of the 


threshold on the properties of the gas or plasma being studied. Interpretation 


in terms of a cascade model will be more difficult because of these factors. 


2 ¢ . e 
For example, the calculations done on this work which seem to agree quit 


well with experimertal results may be good only for the particular laser, 


gas, and optical system used. 


The technique seems to be limited also to low electron densities, 


gas densities, and low plasma temperature. Only further work will tell 


whether these problems can be overcome. 
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It would be interesting to study the breakdown process in the 
presence of weaker discharges such as glow discharges. Here the light emission 


from the plasma may be weak enough so that the problem of seeing the spark 


will be non-existent. The study of the distribution of the photons causing 


ionization emitted by the discharge should be pursued further. 
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APPENDIX 


THIS PROGRAF SOLVES THE RATE FOUATIONS FOR THE RLECTRON DENSITY 
AND THE ELECTRON PPYEERATURE WHICH DESCRIBE TYR BUILDUP IN 

TEE FOCUS OF A LENS DUH TO A LASER BEAM. ENPRGY IS GAINED FROM 
THE RADIATION PLSLD VIA ELECTSON-NEUTRAL INVPRSE SRENSSTRAHLUNG. 
PARTICLE LOSS MECHANISY IS DIPFUSION OP ELECTRONS OUT SO RS rit 
FOCAL VOLUYE. ENERGY IS LOSY DOE TO DIPFUSION, ELASTIC COLLISION 
AND INELASTIC COLLISIONS. RFCOMBINATION,RADIATION, HYDRODYNAMIC 
EPPECTS,COLLISIONS OF THE SECOND KIND, AND IONLZATLON OF EXCITED 
STATES ARE N&SGLECTED TN THE CALCULATION. BREAKDOWN IS DEFINED AS 
THE ATTAINMENT OF AN ELECTRON DENSITY IN THEE FOCUS EQUAL TO 
NB=9.UE13*PK WHERE Ph IS THE GAS PRESSURE IN TORR. 


NT = RLECTRON-ION DENSTYY 

NB = DPNSITY AT BREAKDORN 

DB = BACKGROUND }LFECTKON DENSITY 
DL = DIFFUSION LENGTH 

PR = GAS PRESSURE IN TORR 


EI IONIZATION POTENTIAL IN FV 
EX = EXCITaTION ENERGY IN EV 

PMAX = PEAK LASER POWER IN W/CM*#2 
DIV = CONSTAST 

TPRES = PRESENT TIME 

TT = ELFCTPON THOMPERATUKE IN EV 
ICOUNT = COUNTER 

R = RADIUS UF FOCAL S¥FOT 

A = POCAL ASFA 

ENERGY = LASER PULSE ENERGY 
EXTERNAL GAIN, SINK, ELOSS, HEAT, POWR 
REAL NT, NB 

COMMON DB,DL,PR,EI,EX,PMAX,DIV 
PR=4000. 

PMAX=1.46E8 

NT=1.£E5 

TT=.025 

TPRES=0.0 

DIV=1. 

ICOUNT=1 

NB=9 UE 134#PR 

DL=6 .71-3 

BX=11.5 

EI=15.8 

R=.01 

A=1.57E-4 

ENERGY= 12 102-7*PBAX¥A 

PRINT 101, PMAX, ENERGY 

PRINT1C00, DL,NT,TT,NB,DIV 


NT=NT/DIV 
DB=NT 

NB=N5/D1V Pa 

CALL FD(TT,NT, TPKES, HEAT, ZLOSS, GAIN, SINK, POWR, NB) 

Sone 514.27, F 14 07, EV oT, EM oT, EW 7) 


FORMAT (°0%,5X,E14.7,F 1467, 
FORMAT (80°, *PMAX=",£14 07, ENERGY=" 214-7) 


END 
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SUBROUTINE FD (TT, NT ,TPRES ,WEAT, ELOSS ,GAIN, SINK, POWR, BD) 
SUBROUTINE PD SOLVES Tak SEY OF DIFPEEENTIAL EQUATIONS POR THR 
ELECTRON DENSITY AND TPMPERATURE. TIME STEP SIZE IS FQUAL To W 
CONTROL IS RETUNLNED 10 MAIN PROGRAM IF THE ELECTRON DENSITY 
EQUALS OR UXCERDS NA OF WHEN THE PARTICLE LOSS TERM PXCEEDS 
THE PARTICLE GAIN TENM. FUNCTIONS GXIN (PARTICL!. GAIN), SINK 
(PARTICLi, LOSS), “Sat (..NERGY GAIN), AND ELOSS(EZNERGY LOSS) ARE 
CALLED IN THE CALCULATIONS. 
REAL nT 
COMMON DB,DL,PR,EI,EX,EMAX,DIV 
W=1.E-11 
I=0 
PRINT 100 
TPRES=TPRES4W 
DTEMP=nNT 
TTEMP=TT 
1 POW=PoWk (TPRES, PMAX) 
1=I+1 
DO 10 L=1,2 
D1=GAIN (TTEYP,PR) -SINK (DL, DT2MP,TTEMP, PR) 
CI=HEAT (PK, TTEMP,POW) -FLOSS (TTEMP, DL, FR, EI,EX,DTEMP) -TTEMP*D1 
IF(D1.LT.0.) RETURN 
TTEMP=TL+W*C1 
DTEMP=NT/ (1.-W*D1) 
10 CONTINUE 
NT=DTEMP 
TT=ITEMP 
TPRES=TERES+W 
IF(NT.GT.8D) GO YO 11 
IF (TPRES.GT.3.E-7) RETURN 
IF(I.L1.100) GO TO 1 
I=0 
NT=NT*DIV 
PRINTI01, TPRES,NT,TT,C1,D1 
NT=NT/DIV 
‘GO TO 1 
11 PRINT102, TPRES 
100 FORMA o7°,5X,° / "DENSITY *,5X,"TIMPERATURE*,6X,°EGAIN®, 
‘ (°1°,5X,*TIME®, 9X, *DENSITY °, 5X, 
15X,"PGAIN*) 
101 PORMAT ('0°,5214.7) 
102 FORMAT ('1*,*PREAKDOWN BEFORE THE END OF LASFR PULSE", 5X, 
VeTIME=",£14.7) 
END ee 
SCTION SLOSS (TT,DL,P, LI, 8X,DE 
PUNCTION ELOSS CORPUTES THE ENERGY LOSSES DUE TO IONIZATION, 


EXCITATION, DIPFUSIUN AND ELASTIC COLLISIONS. THERE IS A GRADUAL 
TRANSITION FROM ¥REE DIFFUSION TO AM4IPOLAR DIFFUSION SIMILAR 
TO THE THAT DONE 3Y ALLIS AND ROSE. — 
S = HATIO OP RLECTRON MOBILITY TO ION MOBILITY — 4 
THE RATES HAVE S5EN ObTAINED BY BY AVERAGING TA PROPER EXPRESSIONS 
OVER A MEXWELLIAN DISTRIBUTION. A N FIGHT ORDER POLO. 
YAS BEEN PITTED INTO THESE NUMERICALLY OBTALNED RATES AND 
BFLO® (E1,£2,E3, 54) 
COMMON 0,A,B,C,D,PMAX,DIV 
DOUBLE PRECISION E1,B2,83,E4,T 92¢Se51,22,DL 
DB=Q 
DR=DF*DIV 
7 T=1T 
S=270. 
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K=DE*LDL**«2 7 
TF (LT. WEG) GO PO 30 
IP (X.LT.5.E10) GO TO 4O 
X 1=ALOG (X) 
P =.01944+1./(X1437.26) 
GO TO 50 
30 F =1.0643069# (x** (-.0291024) ) 
GO TO 50 
YO XT=RLOG (X) 
Be en 2119440599980 754% 1- .2350 787 166%X 14X14..008 1446 37934X 1 
F=EX° (P) 
50 CONTINUE 
DE=DE/DIV 
4 IP(T.cT..09) GO To 2 
T EI=-90504.62565+ 18710134 .96*T* (1.-62.0956 13134*T* (1-8. 290271858 
T#T# (104145 .02734454T# (1.-31.971 263874 TH (16-13 2631736347 * (Te 
26 04131925764T* (1.-2.5547973294T))))))) 
EV=41* (1.-DB/DE) 
E£2=159 548574 5-34866. 1747 *T* (1.-70 8520116 8*T* (1.-25 B45H0Z3 RUST ® 
T (1.-10 £6 763389 7*T* (1.41.0220645234*T* (1.-125.96071234T* (1.- 
28.9667 64848 4T* (1.-2..96 36475354%T))))))) 90. 3436225 
13 E3=6 .69U7* ((U.U2F-4FT 411 447 ¥T 41, 058TH 93—. 074 #7) #44) FDEXP (-15.8/T) 
Vt (-4G7.24T-15. VAT HTH. TIS *TEH I+ O74 *THRG) SDEXP (-4O./T) ) / (1#*17.5) 
E4=3.39E 10*DEXKP (1.-11.5/T) /( (1 ##1.5) ®((1.+11.5/T) **2) ) 
GO TO 3 
2 B1=61536 2476154 27525706.664P* (1.—-1.90840 186 3*T* (1.-1.036 365552* 
IT * (1.-—.617 164 29447 * (1.-2 3915656 7281* (1.-.254244 72964T* (1.- 
221644163265 *T* (1.-.1024558964T* (1.-.05797 17 264T* (12-00 250 227374T 
3))))))))) ; 
E1=£1* (1.-D5/DF) 
EB2=4 .97P5* (T#*1.5) #(1./( (1+ 0 045*T) €#3) -16/( (1.4 -115*1) #*3) ) 
1-90.3438225 
16 E3=0.69E7*% ((4o42P-AHTH+ TD GH PETS 1058 THE 3I— 07447 #44) HD EXP (-15 23/7) 
14 (-G07 .2*T-15.18T 47 + 7 35 4T HHS 4074 %T S44) ¥DEXP (-40./T)) / (T**15) 
BU=3.39E 10*DSXE (1 0-1165/T) J ( (T#*# 105) * ( (1241125/T) **2) ) 
3 ELOSS=P*E1/ (PHDL##2) Fi 2* P+ I*E LPs EUS POE y 
ELOSS=ELOSS#2./3. 
RETURN 
END : 
FUNCTION HEAT(P,TT, PW) 
DOUBLF PRECISION G,T 
COMMON DB 
T=TT 
IF (T.GT..09) GO TO 2 
4G = .005262611632-1-476274 564 *T* (12-94 2700 36 7T58*T* (To- 
14 308783698 7#E* (16-25 991445 244 TH (12-15 621324399*T% (1.-9 2201797018 
2*T* (1.-5.1281264 034T* (12-2. 1964 35630*T))))))) 
HEAT=G 
: Pa peel erie. ts) #(1./( (10 4+-G45*T) #43) -16/( (1.4. 1154) ** 3) ) 


HEAT=G 
3 HEAT=HEAT*P*PH*2./3. 
RETURN 
END 
PUNCTION SINK(DL,D%,TT,P) 
COMMON Q,A,B,C,D,PMAX,DIV 
DOUBLF PRECISION A1,2,S,S1,22,T,DL 
DB=Q 
DE=DE*DIV 
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S=270. 

X=D 2*DL¥**2/7 

IPs Lt. 1.86) Co TO 36 

IP(X<LT.5.819) Go To 60 

XT=A LOG (X) 

PF =.019441./ (1437.26) 

GO To 50 

P =1.0693009% (X** (~.0291024) ) 

GO TO 50 

X 1=ALOG (X) 

Pea cr copes 80 Toe Xa 2235078 ANGST ex) + OORINNOS Tolar. 
F=EXP (F) 

CONTINUE 

DE=DE/DIV 

IP (T.GT..09) GO TO 2 

AT=- 440449 .2602493887505.54T* (1.-67.8 196688*T* (12-16 .56 9336 31*T# 
4 (72446 210972895 4T* (1.-38..32996662*T# (1.-14 037828 13*T* (1.- 
26.6022 108 0U4T* (1.-2.600287366*T) )))))) 

A1=A1* (1.-DE/DE) 

GO TO 3 

A1=6071182 £958 44 3257636 .95*T* (1.-2.62388018 941% (1.-1. 1338399834 T* 
1 (1.-.64546 9502*7* (12-240 1604116 4T* (1.—-.258364115"T 4 (1.-. 166187809 
2*T* (1.-—. 103223599 *T* (12-6 0585264 9B8¥T# (12-.025123406*1))))))))) 

A 1=h 1* (1.-DB/D5) 

SINK=A 1*P/ (P* (DL**2) ) 

RETURN 

END 

FONCTION GAIN (TT,P) 

COMMON DB 

DOUBLE FRZCISION G ,T 

T=TT 

G=P*6 .69E7* ((4 4 2E-4 FT 41714 FT OT 41.05 4TH* 3—. 074 THEY) FDEXP (—15.8/ 
IT) + (-4O7 0 2#T-15 6 1H THHE24.739 TH 43+ O74HTHHG) HDEXP(-4O0/T)) / (T** 1-5) 
GAIN=G 
P=0. 

G=P43,3SE TO#DEXP (1.-11.5/T) J ( (1#* 1.5) *( (1.411.5/T) **2) ) 
GAIN=GAIN+F*G 

RETORK 

END : 

PUNCTION POWR(T1,F MAX 

COMMON DB 

T=T1 

IF(T-6.E-8) 1,1,2 

POWR=PMAX¥1/62E-8 

RETURN 

IP(T-1.2F-7) 3,3,4 

POWR=PMAX* (-0.60*1/6.5-8+1.6) 

RETURN 

IF (T-1.65-7) 5,5,6 

POWR=0.4*PPFAX 

RETURN 

POWE=PMAX* (-O.4*T/1.4E-746 2/70) 

RETURN 

END 
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